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SECTION  I 
INTRODUCTION 

One  of  the  major  hazards  to  flight  safety  today  Is  the  Inflight 
Impact  of  birds.  Aircraft  windshields  and  engine  blading  are  especially 
vulnerable  to  damage.  Because  of  the  Importance  of  this  problem,  both 
NASA  and  the  Air  Force  have  extensive  programs  aimed  at  Improving  the 
Impact  tolerance  of  these  critical  components  by  developing  better 
materials  and  better  structural  designs.  However,  this  work  Is  currently 
being  hindered  by  a lack  of  understanding  of  the  actual  Impact  event. 
Fundamental  research  Is  needed  to  define  both  the  loads  generated 
during  Impact  and  the  response  of  various  structures  to  these  loads. 

Work  Is  currently  being  conducted  In  the  Air  Force  Materials 
Laboratory  to  define  bird  Impact  loads.  This  Is  a very  complicated 
problem  since  the  loads  are  actually  coupled  to  the  response  of  the 
structure.  In  order  to  simplify  this,  researchers  have  decldeu  to 
decouple  the  loads  from  the  target  response.  This  Is  accomplished  by 
studying  the  Impact  of  birds  against  a rigid  surface.  The  actual  coupling 
mechanisms  will  be  studies  In  a later  program. 

Another  program  being  conducted  In  the  laboratory  In  parallel  with 
the  bird  loading  program  concerns  the  development  of  a "bird  substitute" 
material.  This  synthetic  bird  would  replace  real  birds  In  actual  Impact 
tests  of  aircraft  components.  Several  materials  are  presently  being 
considered.  Including  gelatin  and  RTV  (room  temperature  vulcanized) 
rubber. 

In  order  to  fulfill  the  requirements  of  both  programs,  It  Is  obvious 
that  a complete  understanding  of  the  fundamental  principles  of  the  Impact 
process  Is  required.  It  seems  reasonable  to  expect  that  the  basic 
mechanisms  of  Impact  of  gelatin  or  RTV  rubber  should  be  similar  to  those 
of  bird  Impact.  This  Is  based  on  the  belief  that  during  impact  each  of 
these  materials  tends  to  behave  as  a fluid.  This  Is  true  for  most  mate- 
rials at  sufficiently  high  velocities.  The  requirements  for  a material 
to  "flow"  Is  that,  the  stresses  generated  during  Impact  substantially 
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exceed  the  strength  of  the  material.  For  many  projectile  materials  such 
as  aluminum  or  steel,  the  stresses  required  to  cause  the  material  to  fail 
are  sufficient  to  cause  local  failure  In  the  target.  For  example,  the 
result  of  a hypervelocity  Impact  Is  a crater. 

However,  there  would  seem  to  be  a large  group  of  materials,  often 
considered  solids,  whose  strength  are  very  low  relative  to  the  strength 
of  many  target  materials.  For  the  purposes  of  this  study,  this  group  of 
materials  will  be  referred  to  as  soft  body  materials.  Soft  body  Impact 
will  be  defined  as  Impact  In  which  the  stresses  generated  substantially 
exceed  the  strength  of  the  projectile  but  are  well  below  the  strength 
of  the  target  material.  This  Implies  that  soft  body  projectiles  will 
flow  upon  Impact  while  the  target  may  see  little  or  no  plastic  deformation. 

Examples  of  soft  body  impact  are  birds  and  tire  fragments  striking  engine 
blading,  hall  striking  aircraft  windshields,  rain  drops  falling  on  a leaf, 
lead  bullets  striking  a steel  plate,  or  even  bugs  hitting  an  automobile 
windshield. 

With  this  In  mind  It  was  felt  that  a general  survey  of  the  state- 
of-the-art  in  soft  body  Impact  would  aid  greatly  in  an  understanding  of 
the  basic  mechanisms  of  bird  impact  as  well  as  the  Impact  of  various 
bird  substitute  materials.  However,  as  seen  in  the  Literature  Survey 
cf  this  report,  the  results  of  the  literature  search  were  limited.  With 
the  exception  of  water  droplet  and  water  jet  impacts,  there  appears  to 
have  been  no  concerted  effort  to  study  the  basic  processes  Involved  In 
soft  body  Impact. 

I 

For  the  past  20  years,  researchers  have  studied  the  Impact  of  ] 

V 

water  droplets  on  & solid  surface  In  order  to  develop  an  understanding  j 

of  the  mechanisms  Involved  In  the  erosion  of  steam  turbine  blades  and  \ 

the  erosion  of  aircraft  and  missile  structures  in  rain.  Much  progress  ! 

has  been  made  In  this  area,  although  there  is  still  controversy  over 

the  exact  amplitude  and  distribution  of  loads,  as  well  as  the  response  j 

of  both  the  liquid  droplet  and  the  Impacted  surface  to  these  loads.  \ 
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Based  on  the  Insights  gained  from  the  literature  search.  It  appears 
possible  to  directly  apily  much  of  the  knowledge  gained  from  fluid  Impact 
studies  to  the  more  general  area  of  soft  body  Impact.  Wf'-.h  this  In  mind, 
the  author  decided  to  undertake  the  present  research  program.  Ihe 
ourpose  of  this  research  was  to  understand  the  basic  processes  Involved 
In  soft  body  Impact.  Both  theoretical  and  experimental  studies  have 
been  conducted  on  the  Impact  of  various  soft  body  materials.  Including 
bird,  gelatin,  and  RTV  rubber. 

A study  of  this  magnitude  should  be  of  great  benefit  since  the 
results  could  be  used  In  many  areas.  As  discussed  above,  the  most 
Important  application  to  the  Air  Force  would  be  In  aiding  In  the  under- 
standing of  bird  Impact  as  well  as  In  the  development  of  a bird 
substitute  material.  However,  there  are  many  other  practical  applications. 
For  example,  knowledge  from  this  program  would  aid  In  the  study  of  hall 
damage,  rain  erosion,  and  damage  due  to  debris  Impact  during  a hurricane 
or  tornado. 
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SECTION  II 
LITERATURE  SURVEY 

Th.2  impact  of  low  strength  materials  has  only  been  studied  exten- 
sively In  the  last  few  years.  Thus,  there  has  been  very  little  research 
in  the  area  of  soft  body  impact  and  almost  none  concentrated  on  the 
actual  mechanisms  involved  in  the  impact  process.  Early  researchers 
such  as  Hodgkinson  (Reference  1),  Tait  (Reference  2),  Vincent  (Reference 
3),  and  Raman  (Reference  4)  studied  the  coefficient  of  restitution 
(ratio  of  the  rebound  velocity  to  the  initial  velocity)  of  various 
soft  objects.  They  found  that  the  coefficient  of  restitution  decreased 
with  Increasing  velocity,  demonstrating  that  the  percentage  of  kinetic 
energy  dissipated  during  Impact  Increased  with  increasing  velocity. 

Very  little  work  has  been  done  in  this  area  except  for  recent  studies  of 
birds  and  bird  substitute  materials. 

In  the  1960's,  several  English  engine  people,  including  Allcock 
and  Collin  (Reference  5),  studied  the  deflection  of  beams  due  to  the 
impact  of  various  soft  objects  including  birds.  Based  on  a comparison 
of  the  results,  they  determined  that  gelatin  projectiles  and  birds 
produced  comparable  loads  during  impact.  Tudor  (Reference  6),  studying 
the  impact  of  gelatin  projectiles  on  a cantilevered  beam,  developed  a 
functional  relationship  between  the  beam  deflection  and  the  Initial 
momentum  of  the  projectile.  In  America  work  was  carried  out  by  Tsai, 
et  al.  (Reference  7)  who  found  a similar  relationship  for  the  impact  of 
RTV  rubber  balls.  The  only  theoretical  study  of  bird  loading  was 
performed  by  MacCauley  (Reference  8)  and  Mitchell  (Reference  9)  in 
Canada.  The  work  of  both  men  were  only  crude  approximations  of  the  real 
impact  case.  MacCauley  assumed  the  bird  to  behave  as  a perfect  fluid 
and  Mitchell  assumed  the  bird  to  behave  as  a semi-rigid  projectile. 

With  other  simplifying  assumptions,  both  men  derived  expressions  for  the 
pressures  generated  during  impact.  No  experimental  work  was  performed 
to  check  their  results. 
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Since  so  little  has  been  done  In  the  specific  area  of  soft  body 
Impact,  It  was  felt  that  a survey  of  the  general  area  of  Impact  might 
be  of  benefit.  Specific  Interest  was  paid  to  studies  concentrating  on 
the  response  of  the  projectile  to  the  Impact,  as  opposed  to  the  response 
of  the  target.  It  was  found  that  numerous  studies  of  the  mechanics  of 
the  Impact  event  and  the  related  Impact  strength  of  materials  have  been 
conducted  over  a large  range  of  Impact  velocities.  Much  of  this  work  Is 
summarized  In  textbooks  by  Goldsmith  (Reference  10),  Kolsky  (Reference  11), 
and  Johnson  (Reference  12). 

For  the  Impact  of  solids,  the  mechanisms  dominating  the  Impact 
process  vary  with  Impact  velocity.  Realizing  this,  several  Individuals 
have  attempted  to  divide  the  field  of  Impact  of  solids  Into  various 
physical  regimes.  Hopkins  and  Kolsky  (Reference  13)  defined  five  regimes; 

(a)  elastic  Impact,  (b)  plastic  Impact,  (c)  hydrodynamic  Impact,  (d)  Impact 
at  sonic  velocities,  and  (e)  explosive  Impact.  In  elastic  Impact,  the 
stresses  generated  do  not  exceed  the  yield  strength  of  the  material. 

Thus,  the  nature  and  duration  of  the  Impact  only  depend  on  the  elastic 
moduli  and  the  elastic  wave  velocities  of  the  material.  With  Increased 
velocity,  the  stresses  generated  cause  plastic  deformation.  The  material 
strength  Is  still  a dominating  factor.  As  velocities  Increase  still 
further,  the  stresses  generated  by  deceleration  of  the  projectile  greatly 
exceed  the  yield  stress.  In  this  hydrodynamic  regime,  the  projectile  and 
target  can  be  treated  as  fluids,  and  It  Is  the  material  density  which 
dominates  Instead  of  strength.  In  each  of  these  regimes,  stress  waves  I 

propagate  Into  the  material,  thus  dissipating  energy  throughout  the  I 

projectile  and  target.  As  the  velocity  of  Impact  approaches  or  exceeds 
the  wave  velocity,  more  energy  must  be  dissipated  In  the  local  region  of 

Impact.  Wave  motion  plays  an  Increasingly  Important  role  In  determining  ■ 

the  local  stress  distribution.  Shock  waves  are  also  generated.  As  the  ] 

velocity  continues  to  Increase,  all  of  the  energy  Is  deposited  In  the 

local  area.  The  heat  produced  Is  concentrated  In  a very  small  region,  and  : 

Is  thus  sufficient  to  melt  and  eventually  to  vaporize  the  material. 

This  process  Is  analogous  to  a small  explosion  taking  place  on  the  target 
surface. 
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Although  these  regimes  have  been  defined  for  the  Impact  of  high 
strength  projectiles,  there  Is  no  reason  to  believe  that  these  same 
regimes  would  not  be  applicable  to  the  Impact  of  low  strength  materials. 
Barber  (Reference  14)  found  that  cylindrical  projectiles  of  RTV  rubber 
and  gelatin  bounced  at  low  velocities  and  "flowed"  at  high  velocities. 
Similar  results  were  also  found  for  real  birds.  Thus,  a further  look 
into  these  various  regimes  should  lead  to  a better  understanding  of  soft 
body  Impact. 

Johnson  (Reference  12)  described  initial  studies  of  elastic  Impact 
which  applied  elementary,  one-dimensional,  elastic  wave  theory  to  the 
Impact  of  cylinders.  Using  this  simple  theory  to  study  the  problem  of  a 
cylinder  impacting  a rigid  plate,  an  interfacial  pressure  P <*  pcAuQ  and  a 
duration  T - 2L/uQ  can  be  derived,  where  c£  is  the  longitudinal,  elastic 
wave  speed  In  the  projectile.  However,  this  simple  theory  neglects  the 
effects  of  transverse  strain  and  Inertia  as  well  as  all  dissipative 
forces  In  the  projectile  material. 

Theories  for  elastic  Impact  which  take  into  account  the  effects  of 
radial  Inertia  were  first  proposed  by  Pochhamer  (Reference  15)  and 
Chree  (Reference  16)  independently.  Later,  theories  were  proposed  by 
Love  (Reference  17)  and  Rayleigh  (Reference  18)  and  more  recently  by 
Redwood  (Reference  19).  Using  Love's  theory,  Conway  and  Jakubowski 
(Reference  20)  analyzed  the  coaxial  impact  of  nearly  perfect,  square- 
ended,  finite  length  bars.  Their  results  showed  a slight  change  in  pulse 
shape  due  to  the  transverse  effects. 

Donnell  (Reference  21)  In  1930  was  apparently  the  first  to  investi- 
gate longitudinal  plastic  wave  propagation.  However,  major  contributions 
to  the  theory  of  plastic  wave  propagation  were  not  made  known  until 
after  World  War  II  when  Taylor  (Reference  22),  Rakhmatulin  (Reference  23), 
and  Karman  and  Duwez  (Reference  24)  published  their  theories  which  had 
been  developed  independently  during  the  war.  In  a later  paper,  Taylor 
(Reference  25),  in  an  attempt  to  develop  a simple  method  for  determining 
the  dynamic  yield  stress  of  materials,  considered  the  propagation  of  both 
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elastic  and  plastic  waves  during  Impact  of  a cylinder  on  a rigid  wall. 
His  simple,  theoretical  model  led  at  least  to  qualitative.  If  not 
quantitative,  understanding  of  projectile  mushrooming  upon  Impact. 
Experiments  by  Whlffen  (Reference  26)  showed  that  Taylor's  method  was 
Inaccurate  except  over  a limited  range  of  velocities.  Hawkyard  (Ref- 
erence 27)  examined  the  consequences  of  establishing  an  energy  balance 
across  the  discontinuity  at  the  plastic  wavefront.  His  predictions  of 
the  f*nal  shape  of  the  mush-oomed  projectile  corresponded  much  more 
nearly  to  the  shapes  found  experimentally  by  Whlffen  than  did  the  pre- 
dictions of  Taylor. 

Although  the  theories  of  Taylor  and  Hawkyard  aided  in  understanding 
the  mechanisms  of  deformation  during  Impact,  they  were  too  Idealized  to 
give  good  quantitative  results.  Both  effectively  considered  the  pro- 
jectile to  be  a rigid,  perfectly  plastic  material.  In  an  attempt  to 
overcome  this  limitation,  Barenblatt  and  Ishllnskil  (Reference  28) 
considered  the  Impact  of  a linear,  viscoplastic  rod  on  a rigid  wall. 

Ting  (Reference  29)  extended  their  work  to  study  the  Impact  of  a non- 
linear, viscoplastic  projectile  and  contrasted  his  results  with  those  of 
Barenblatt  and  Ishllnskil.  All  of  this  work  considered  one-dimensional 
wave  propagation. 

When  the  stresses  generated  during  impact  greatly  exceed  the  yield 
stress  of  the  material,  the  problem  can  be  approached  hydrodynamlcally. 
In  the  hydrodynamic  approach  used  to  study  high  velocity  Impacts,  the 
material  strength  and  viscosity  are  neglected  and  a simple  pressure- 
density-energy  equation  of  state  Is  used  to  describe  the  material.  For 
most  materials  this  approach  Is  applicable  over  a wide  stress  range. 

Hopkinson  (Reference  30),  In  studying  the  Impact  of  a lead  bullet 
on  a steel  plate,  found  experimentally  that  the  stress  in  the  target  was 
approximately  that  to  be  expected  from  a fluid  jet  whose  density  was  the 
same  as  for  lead.  Blrchoff,  et  al.  (Reference  31)  In  their  study  of 
explosives  with  lined  cavities,  may  have  been  the  first  to  Incorporate 
into  a mathematical  theory  the  Idea  that  a high  strength  material  such 
as  steel  could  be  treated  as  though  It  were  a perfect  fluid  when  sub- 
jected to  very  high  stress  levels. 
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Thus,  It  is  seen  that  even  high- strength  materials  "flow"  when 
subjected  to  very  high  pressures.  The  major  drawback  to  this  area  of 
research  is  that  the  stresses  generated  at  the  projectile-target  inter- 
face are  so  great  that  the  target  material  also  behaves  hydrodynan.ically, 
resulting  in  large  deformation^  (craters)  in  the  target.  Thus,  it  is 
not  possible  to  uncouple  the  response  of  the  projectile  from  that  of 
the  target.  A review  of  the  literature,  as  performed  by  Chakrapani  and 
Rand  (Reference  32),  for  example,  showed  that  most  of  the  research  in 
this  regime  is  concerned  with  predicting  or  measuring  crater  shape  and 
depth.  Therefore,  it  was  decided  that  further  study  of  this  regime  or 
the  higher  velocity  regimes  (sonic  and  explosive  regimes)  would  be 
fruitless. 

A much  more  promising  area  of  research  appears  to  be  the  study  of 
the  Impact  of  fluids,  in  which  the  fluid  projectile  is  always  assumed 
to  behave  hydrodynamical ly.  An  extensive  amount  of  theoretical  and 
experimental  work  has  been  conducted  on  the  impact  of  fluids  on  solid 
targets,  in  vhich  the  target  Is  assumed  to  have  little  or  no  plastic 
deforma .ion. 

The  majority  of  work  In  this  area  has  been  focused  on  the  impact  of 
water  jets  and  water  drop’ets,  both  In  the  study  of  rain  and  steam 
erosion  as  well  as  water  jet  cutting.  In  these  Investigations,  there  are 
several  features  of  tho  impact  process  of  greatest  importance.  Heymann 
(Reference  33)  stated  that  these  were;  (a)  the  initial  impact  pressure, 

(b)  the  area  over  which  the  pressure  acts,  (c)  the  velocity  of  the  lateral 
outflow  of  liquid  after  impact,  and  (d)  the  approximate  decay  time  of 
high  Impact  forces. 

Some  of  the  early  work  was  conducted  by  Cook  (Reference  34)  in  1928, 
in  which  he  compared  the  impact  of  a column  of  water  to  the  waterhammer 
problem.  He  stated  that  in  both  cases,  a moving  column  of  water  is 
suddenly  an^sted  by  a fixed  surface,  thereby  causing  an  Instantaneous 
pressure  In  the  front  layer  of  the  column  given  by  the  expression 
P - uQ  v£7T,  where  uQ  is  the  Impact  velocity,  p is  the  initial  density 
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of  the  water,  and  8 Is  Its  compressibility.  Noting  the  magnitude  of  this 
Initial  "waterhammer"  pressure,  he  stated  that  "the  pressi'r^ 
on  an  “le.Tent  of  :u. 'acc  **  its  first  encounter  with  water  moving  at  a 
finite  velocity  Is  different  from  that  produced  by  the  steady  Impact 
of  a moving  stream  of  water  at  the  same  velocity".  Heymann  (Reference  35) 
phrased  It  differently,  noting  that  In  this  type  of  Impact,  compressibility 
effects  rather  than  flow  phenomena  Initially  predominate  In  the  liquid. 
Thus,  Cook  and  Heymann  Identified  the  two  basic  regimes  of  water  jet 
Impact,  the  Initial  shock  regime  and  the  late-time  steady  flow  regime. 

These  two  regimes  were  later  Identified  by  Bowden  and  Brunton 
(Reference  36)  In  their  more  detailed  description  of  the  impact,  of  a 
cylinder  against  a rigid  plate.  In  their  paper,  they  stated  that  upon 
Impact  of  an  Ideal,  flat-ended  cylinder  on  a rigid  plate,  stress  waves 
propagate  Into  the  liquid  and  the  pressure  at  the  Interface  Is  the 
waterhammer  pressure  given  by  Cook  (Reference  34),  P « pcQu0,  described 
In  terms  of  the  wave  speed,  cQ,  Instead  of  the  compressibility, 3 . 

According  to  their  description,  as  release  waves  propagate  Inward  to  the 
center  of  the  cylinder,  liquid  Is  turned  outward  and  flows  radially  with 
the  same  velocity  uQ  and  the  pressure  at  the  Interface  Is  released  and 
decays  to  the  steady  flow  pressure.  The  pressure  reduces  to  zero  when 
all  of  the  liquid  has  been  turned.  The  duration  of  the  high  pressure 
was  determined  by  the  time  taken  for  the  release  waves  from  the  sides  of 
the  liquid  cylinder,  traveling  at  velocity  cQ,  to  reach  the  center  of  the 
cylinder. 

For  the  Impact  of  a sphere  on  a rigid  plate,  Bowdes,  and  Field 
(Reference  37)  described  the  event  to  be  somewhat  the  same.  Upon  Initial 
impact,  the  pressure  at  the  Interface  again  reaches  a maximum  given  by 
the  waterhammer  equation.  Release  waves  propagate  inward  sooner  than 
for  a cylinder  of  equal  radius  so  that  the  duration  of  the  peak  pressure 
Is  less.  They  also  stated  that  the  waterhammer  equation  was  only  correct 
for  low  velocity  Impacts.  For  high  values  of  impact  velocity,  the  pres- 
sures are  such  that  cQ  Is  more  correctly  replaced  by  the  shockwave 
velocity,  us.  Thus,  they  stated  that  the  Initial  pressure  generated 
upon  impact  of  a cylinder  or  sphere  is  given  by  the  relationship 
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■ pusu0.  (This  pressure  will  be  referred  to  as  the  Hugonlot  pressure, 
since  It  can  be  derived  from  the  Hugonlot  relationships  across  a shock.) 

Although  these  two  models  of  cylinder  and  sphere  Impact  are 
generally  accepted,  there  Is  much  disagreement  In  the  literature  over 
certain  features  of  the  Impact  event.  Host  of  this  disagreement  centers 
around  the  peak  pressures  generated  Initially  during  Impact  and  the 
Initial  radial  velocity  of  the  released  liquid. 

It  Is  generally  agreed  that  Bowden  and  Field  (Reference  37)  were 
correct  In  stating  that  the  Instantaneous  pressure  generated  during 
impact  of  a cylinder  is  the  Hugonlot  pressure.  However,  oreat  difficulty 
has  arisen  In  verifying  this  experimentally  due  to  the  inability  to 
generate  a flat-ended  cylinder  of  water  which  can  be  impacted  on  a rigid 
surface.  Bowden  and  Brunton  (Reference  36)  described  a technique  for 
accelerating  a jet  of  water  against  a plate.  However,  these  jets  suffer 
from  an  Inherently  rounded  front  surface.  Both  Brunton  (Reference  38) 
and  Kinslow,  et  al.  (Reference  39)  were  unable  to  obtain  Instantaneous 
peak  pressures  upon  impact,  probably  a result  of  the  rounded  front. 

Thus,  Kinslow  recorded  a peal,  pressure  of  only  half  that  predicted. 

Field,  et  al . (Reference  40),  also  recorded  peak  pressures  less  than 
the  Hugonlot  pressure.  Thus,  work  continues  to  develop  better  methods 
of  testing  for  cylindrical  Impact. 

Unlike  the  case  of  cylinder  Impacts,  there  's  much  less  agreement 
concerning  the  magnitude  and  distribution  of  the  initial  pressure  in 
spherical  droplet  impacts.  Engel  (Reference  41),  an  early  pioneer  In  the 
field  of  erosion,  stated  that  the  average  pressure  in  the  region  of  the 
droplet  which  had  been  traversed  by  the  shock  could  be  expressed  by 
P - | Pc0uo  where  Q*  * constant,  approaches  unity  for  high  velocity 
Impact  and  the  1/2  Is  a consequence  cf  the  spherical  shape  of  the  droplet. 
However,  she  made  no  attempt  to  specify  Instantaneous  distribution  of 
pressure.  Jenkins  and  Booker  (Reference  42)  agreed  with  Bowden  and 
Field  (Reference  37)  that  the  maximum  Impact  pressure  for  a curved  surface 
Is  the  Hugonlot  pressure,  and  that  this  pressure  Is  uniform  over  the 
liquid-solid  interface  until  lateral  jetting  begins.  In  contrast  to 
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this,  Skalak  and  Felt  (Reference  43)  concluded  that  the  "average"  pressure 
over  the  llqu’d-solld  Interface  was  the  Hugonlot  pressure  but  also  con- 
cluded that  the  pressure  was  no',  uniform,  with  maximum  pressure  (greater 
than  pusuq)  near  the  edge  of  contact.  Heymann  (Reference  44)  developed 
a theory  based  on  the  assumption  that  oblique  shocks  exist  at  the  droplet 
edge  and  used  this  theory  to  predict  that  the  pressure  Increases  during 
early  Impact.  This  Increased  pressure  exists  at  the  droplet  edge  and 
continues  to  Increase  until  lateral  jetting  starts.  Pressures  exceeding 
twice  the  Hugonlot  pressure  were  predicted  for  an  Impact  of  Mach  • 0.3. 

Even  greater  disagreement  resulted  frum  contradicting  experiments. 

As  stated  earlier,  the  jet  launching  technique  used  by  Brunton  (Ref- 
erence 38)  produced  a cylinder  with  a rounded  Impact  face.  Thus,  this 
jet  was  used  by  several  researchers  to  simulate  droplet  Impact.  However, 
Brunton's  experiment  permitted  the  determination  of  average  pressure 
only,  Instead  of  the  actual  distribution.  Using  this  same  technique, 
Klnslow,  et  al . (Reference  39),  found  the  maximum  pressure,  about  half 
the  predicted  Hugonlot  pressure,  to  occur  at  the  center  of  the  Impact. 
Johnson  and  Vickers  (Reference  45)  found  approximately  the  same  pressure 
at  the  center  but,  contrary  to  the  work  of  Klnslow,  et  al.,  measured 
maximum  pressures  at  the  edge  greater  than  the  Hugoniot  pressure.  Using 
a two-dimensional  droplet  impact  test,  Rochester  and  Brunton  (Reference 
46)  verified  the  work  of  Johnson  and  Vickers.  Using  high  speed  "shadow- 
graph" photography  which  allowed  them  to  observe  the  stress  waves 
generated  In  the  droplet  during  Impact,  Rochester  and  Brunton  observed 
a region  of  high  pressure  located  at  the  edge  of  drop  prior  to  lateral 
jetting.  Field,  et  al.  (Reference  40),  used  an  experimental  technique 
of  launching  a target  projectile  at  a suspended  droplet.  Using  a 
similar  high  speed  photographic  system  to  that  of  Rochester  and  Brunton, 
they  also  observed  a region  of  high  pressures  near  the  edge  of  the 
droplet. 

Not  only  Is  the»*e  disagreement  In  the  literature  over  the  amplitude 
and  distribution  of  peak  pressures  generates  during  initial  stages  of 
Impact,  but  there  is  also  disagreement  over  the  lateral  flow  velocities 
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resulting  during  these  Initial  stages.  Using  the  various  experimental 
techniques,  Engel  (Reference  41),  Jenkins  and  Booker  (Reference  42), 

Fyall  (Reference  47),  Camus  (Reference  48),  and  Rochester  and  Brunton 
(Reference  46)  observed  radial  velocities  from  three  to  seven  times  the 
Initial  Impact  velocity  of  the  jet  or  droplet.  Bowden  and  Brunton 
(Reference  36)  postulated  that  the  high  radial  velocity  resulted  from 
jetting  action  of  the  water  between  the  confines  of  the  drop  and  the 
Impacted  surface.  The  results  of  Rochester  and  Brunton  (Reference  46) 
correlated  well  with  this  theory.  Contrary  to  this,  Jenkins  and  Booker 
(Reference  42)  found  that  their  results  compared  favorably  with  the 
theory  that  the  radial  wash  originates  as  steady,  incompressible 
Bernoulli  flow,  where  the  maximum  stagnation  pressure  1r  the  initial 
peak  Impact  pressure.  The  experimental  work  of  Fyall  (Reference  49) 
agreed  with  this  theory.  Even  studies  of  ice  Impact  by  McNauqhton, 
et  al.  (Reference  50)  showed  this  theory  to  accurately  predict  Initial 
flow  of  the  Ice  particles  upon  impact. 

Besides  development  and  refinement  of  experimental  techniques,  a 
large  effort  has  been  applied  to  the  use  of  numerical  techniques  to 
develop  a better  understanding  of  the  impact  event.  Huang,  et  al. 

(Reference  51)  spent  a considerable  amount  of  time  in  the  study  of 

cylindrical  and  spherical  impact  at  subsonic  velocities.  For  both  of  ? 

these  cases,  they  found  the  maximum  pressure  to  occur  at  the  center  of 

impact  and  to  have  a magnitude  less  than  the  Hugoniot  pressure.  However, 

these  results  are  in  doubt  since  the  numerical  technique  was  unable  to 

properly  handle  the  strong  shocks  generated.  The  apparent  introduction 

of  artificial  viscosity  probably  prevented  the  program  from  responding  1 

rapidly  enough  to  measure  the  actual  peak  pressures.  Also,  the  numerical  j 

solution  predicted  precursor  waves  upon  impact  which  propagated  at 
velocities  much  greater  than  the  shock  velocity.  However,  there  is  no 
apparent  physical  reason  for  the  existence  of  these  waves. 

Glenn  (Reference  52)  did  similar  studies  in  order  to  point  up  some 
of  the  discrepancies  of  the  work  by  Huang,  et  al.,  and  to  extend  the 
scope  to  include  Initially  supersonic  impact.  All  of  his  work  concentrated 
on  the  Impact  of  a flat  cylinder  on  a rigid  plate.  For  both  subsonic  and 
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supersonic  impact,  the  peak  pressure  at  the  surface  was  the  Hugonlot 
pressure,  thereby  throwing  doubt  on  the  work  by  Huang,  et  al.  Glenn 
also  found  the  maximum  radial  velocities  to  be  greater  than  the  Initial 
Impact  velocity.  Thus,  Glenn's  results  compared  well  with  experimental 
results.  Numerical  studies  by  Lesser  and  Field  (Reference  53)  and 
Prichett  and  Riney  (Reference  54)  also  predicted  the  Initial  peak  pressure 
at  the  center  of  a cylinder  to  be  the  Hugonlot  pressure. 

As  stated  earlier,  once  radial  release  waves  have  propagated  Into 

the  center  of  the  cylinder  (or  sphere)  the  pressure  begins  to  decrease. 

After  several  wave  reflections,  the  flow  of  a cylinder  of  fluid  will 

approach  a steady  state  condition.  At  the  center,  the  pressure  will  be 

equal  to  the  stagnation  pressure  of  a jet.  From  Bernoulli's  Equation, 

2 

the  stagnation  pressure  Is  seen  to  be  P = l/2puQ  for  an  Incompressible 
fluid.  For  a low  velocity  impact,  this  pressure  is  an  order  of  magni- 
tude less  than  the  maximum  shock  pressure  occurring  initially  at  the 
center. 

The  major  area  of  study  in  the  steady  flow  regime  is  the  distribution 
of  pressure  and  velocity  over  the  target  surface.  The  pressure  and 
velocity  are  related  by  Bernoulli's  Equation.  Thus,  knowledge  of  one 
implies  knowledge  of  the  other.  However,  it  is  not  possible  to  calcu- 
late the  distribution  of  pressure  (or  velocity)  in  the  impact  region  for 
an  axi symmetric  jet.  An  exact  solution  can  be  obtained  for  a two- 
dimensional  jet.  Milne-Thompson  (Reference  55)  derived  the  following 
expression  for  the  potential  flow  velocity  distribution  along  the 
surface  of  a plane  jet  impinging  on  a flat,  rigid  surface: 
y/b  = 2/tt  [arctan  (u/uQ)  + arctanh  (u/uQ)]  where  y is  the  distance  from 
’ the  center  along  the  surface,  b is  the  initial  width  of  the  jet,  u is 
the  velocity  along  the  surface,  and  uQ  is  the  initial  impact  velocity. 

This  solution  assumes  a free-streamline  jet,  thereby  implying  that  at 
large  distances  from  the  center  of  impact,  u = u_.  This  expression  is 
substituted  into  Bernoulli's  Equation  of  the  form  P - PQ  = 1/2  puQ 
[1  - (u/uQ)2].  Using  elementary  momentum  considerations.  Banks  and 
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Chandrasekhara  (Reference  56)  developed  a normal  curve  approximation  to 
the  exact  pressure  distribution  of  the  form  P - P0  * 1/2  puQ2  exp  [jj-  (j£)]. 
This  distribution  closely  resembles  that  of  Mllne-Thompson.  Using 
a similar  approach.  Banks  and  Chandrasekhara  developed  ar.  error  curve 
approximation  for  an  axisymmetrlc  jet  Impact  of  the  form 
P - P » 1/2  pu*  exp  ( - (7-) 2]  where  r Is  the  radial  distance  from 

0 0 9 

the  center  and  a Is  the  Initial  radius  of  the  jet.  This  expression  was 
a fair  approximation  to  the  data  obtained  by  Gibson  (Reference  57), 

Leach  and  Walker  (Reference  58)  found  that  their  experimental  data 

? 

was  fit  well  by  a polynomial  expression  of  the  fo^'m  P - P *■  1/2  puQ 
(1 -3(^)2  + 2(^)3)  where  R was  the  value  of  r where  P - PQ.  They 
found  from  momentum  considerations  that  R/a  ~ 2.58.  From  this  It  was 
concluded  that  the  region  where  the  pressure  Is  significantly  greater 
than  the  ambient  pressure  is  confined  to  about  2.6  jet  radii. 

Taylor  (Reference  59),  studying  oblique  impact,  stated  that  the 
total  downward  force  on  the  target  plane  decreases  as  the  angle  of 

Incidence  decreases.  However,  the  maximum  (stagnation)  pressure  is  the 

2 

same  for  all  angles  of  impact,  namely  1/2  pu  , although  the  area  over 
which  this  maximum  pressure  acts  Is  much  reduced  as  the  angle  of 
Incidence  becomes  small. 

Based  on  the  above  information,  several  points  are  clear.  There 
is  general  agreement  about  the  pressures  generated  during  the  impact 
of  a cylinder  of  water  on  a rigid  plate.  The  pressure  at  the  center  is 
initially  pusuQ,  the  Hugonlot  pressure,  and  remains  at  the  level  until 
lateral  release  waves  permit  radial  flow  to  begin.  After  a sufficient 
time,  steady  flow  Is  set  up  with  the  pressure  at  the  center  given  by 
Bernoulli's  Equation. 

There  Is  less  agreement  about  the  Impact  of  a sphere  (or  droplet) 
of  water.  Physically,  one  expects  the  initial  pressure  at  the  center 
to  be  the  same  as  for  a cylinder,  pu$u0,  and  to  remain  at  this  level 
until  radial  release.  In  addition  there  is  extensive  experimental 
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evidence  that  as  the  droplet  Impinges  on  the  target,  the  peak  pressure 
moves  outward  with  the  point  of  contact  and  Increases  In  magnitude  until 
lateral  jetting  (radial  flow)  begins.  The  duration  of  the  droplet  Impact 
Is  probably  too  short  to  allow  sufficient  time  for  steady  flow  to  be 
completely  established. 
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SECTION  III 
IMPACT  THEORY 

In  the  following  section,  hydrodynamic  theory  Is  applied  to  the  \ 

impact  of  a soft  body  against  a flat  target.  Equations  are  derived  for 
the  pressures  generated  during  Impact.  Much  of  this  theory  Is  based  on 
the  Ideas  presented  In  the  Literature  Survey.  Concepts  from  various 

* 

disciplines  have  been  combined  to  develop  a comprehensive  Impact  theory. 

Initially,  the  normal  Impact  of  a right  circular  cylinder  onto  a rigid 
plate  is  considered.  Next,  several  special  cases  are  considered. 

Including  oblique  impact.  Impact  of  a yawed  projectile.  Impact  of  a 
projectile  with  a curved  leading  edge,  addition  of  porosity  to  the 
projectile  material,  and  the  Impact  onto  a non-rigid  target.  A subsection 
is  included  which  explains  the  development  of  the  equations  of  state 
used  to  represent  soft  body  materials.  A significant  portion  of  this 

work  deals  with  altering  the  equations  of  state  to  account  for  the  j 

presence  of  porosity.  Finally,  the  limitations  of  this  impact  theory  are 

presented. 

" 

1.  HYDRODYNAMIC  THEORY  j 

| 

When  a projectile  of  any  material  Impacts  a target  plate,  the  j 

particles  on  the  front  surface  of  the  projectile  are  instantaneously 

brought  to  rest  relative  to  the  target  face  and  a shock  is  formed.  The  1 

purpose  of  this  shock  wave  is  to  bring  each  succeeding  layer  of  particles  j 

to  rest.  The  shock  compression  of  a layer  of  particles  is  so  rapid  that  \ 

the  particles  away  from  the  edge  of  the  projectile  do  not  have  time  to  ] 

I 

"communicate"  with  the  free  surface.  This  implies  that  these  particles  } 

behave  as  If  they  are  In  a semi-infinite  medium  which  can  undergo  only 
plane  strain  compression.  Thus,  shock  compression  In  a bounded  medium 
Is  usually  considered  to  be  a plane  strain  process. 

The  pressure  In  the  shock  compressed  region  Is  very  high  initially 
and  Is  constant  throughout  the  region  at  early  times.  As  the  shock 
propagates  up  the  projectile,  the  particles  along  the  projectile's  edge 
are  subject  to  a very  high  pressure  gradient  due  to  the  shock  loading 
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on  one  side  and  the  free  surface  on  the  other.  This  pressure  gradient 
causes  the  particles  to  be  accelerated  radially  outward  and  a release 
wave  Is  formed.  The  function  of  this  release  wave  Is  to  relieve  the 
radial  pressures  In  the  projectile. 

A very  complicated  state  of  stress  begins  to  develop  In  the  pro- 
jectile. The  radial  pressure  release  causes  shear  stresses  to  appear. 

The  radial  acceleration  of  the  particles  also  results  in  tensile 
stresses  being  developed.  If,  at  any  time,  the  state  of  stress  Is  such 
that  the  strength  of  the  material  Is  exceeded,  the  material  will  "flow". 

For  soft  body  Impact,  It  will  be  assumed  that  the  stresses  throughout 
the  Impact  event  greatly  exceed  the  material  strength,  so  that  the  flow 
will  continue  Indefinitely.  For  these  materials,  to  a first  approximation 
the  material  strength  can  be  neglected  so  that  they  can  be  considered  to 
behave  as  fluids. 

After  several  reflections  of  the  release  waves,  a condition  of 
steady  flow  Is  established.  A constant  pressure  and  velocity  field  Is 
set  up  In  the  projectile,  and  the  particles  flow  along  paths  which  are 
fixed  In  space,  called  streamlines. 

In  the  following  section,  a hydrodynamic  impact  theory  will  be 
used  to  study  the  mechanisms  Involved  In  soft  body  impact.  In  order  to 
simplify  the  analysis,  the  case  of  e homogeneous,  right-circular 
cylinder  impacting  normally  on  a rigid  plate  will  be  considered  initially. 
By  assuming  the  target  to  be  rigid,  the  response  of  the  target  material 
to  the  Impact  pressures  can  be  neglected.  This  simplifies  the  problem 
and  allows  the  researcher  to  concentrate  on  the  response  of  the  projectile. 
Several  assumptions  are  implicit  In  this  theory.  The  projectile  material 
will  be  considered  homogeneous,  even  when  high  amounts  of  porosity  are 
considered.  As  stated  previously,  the  strength  of  the  projectile  material 
will  be  neglected.  Also,  In  order  to  simplify  the  analysis,  the  material 
viscosity  and  the  frictional  forces  along  the  impact  surface  will  be 
neglected. 
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a.  Shock  Regime 

For  the  normal  Impact  of  a cylinder  on  a rigid  pUte,  the  flow 
across  a shock  can  be  considered  one-dimensional  and  adiabatic,  Irrevers- 
ible. Figure  1-b  Illustrates  a shock  wave  propagating  Into  a fluid  at 
rest,  where  u$  Is  defined  as  the  velocity  of  the  shock  propagating  Into 
the  fluid  at  rest  and  up  Is  the  velocity  of  the  particles  behind  the  shock 
In  this  reference  system.  From  this  figure,  it  can  be  seen  that  the 
particle  velocity  Is  actually  the  change  In  velocity  across  the  shock. 

Figure  1-b  Illustrates  the  case  for  which  the  velocities  are  all 
measured  relative  to  the  fluid  In  the  shocked  state.  This  case  is  syn- 
onymous with  the  Impact  of  a cylinder  on  a rigid  plate.  The  projectile's 
Initial  velocity  Is  seen  to  be  uQ  and  it  Is  brought  to  rest  behind  the  shock. 

In  order  to  write  the  conservation  laws  across  the  shock,  the  steady 
state  shock  condition  must  be  considered.  This  case  Is  shown  In  Figure  1-c. 
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Figure  1.  One  Dimensional  Shock  Flow,  a)  Shock  Propagating  Into 
a Fluid  at  Rest,  b)  Flow  Brought  to  a Rest  Across 
the  Shock,  c)  Standing  Shock. 
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For  this  case*  the  equations  of  conservation  of  mass  (continuity) 
and  momentum  may  be  written 

pius  = P2(VV  (1) 

pi  + pius2  = p2  ♦ p2(us-up)2  <2> 

Combining  these  two  equations,  the  pressure  behind  the  shock  Is  found 
to  be 

Vpi = piVP  (3) 

The  pressure  In  the  shocked  region,  given  by  Equation  3,  Is  often 
referred  to  as  the  Hugonlot  pressure.  Throughout  the  remainder  of 
this  report,  this  pressure  will  be  represented  by  the  notation,  PH. 

For  the  impact  of  a cylinder  on  a rigid  plate.  It  can  be  seen  that 
Up  ■ uQ.  Thus,  for  this  case.  Equation  3 becomes 

PH  8 piV0  00 

For  very  low  impact  velocities,  the  shock  velocity,  u , can  be 
approximated  by  the  Isentroplc  wave  speed  in  the  material,  cQ.  Thus, 
for  low  impact  velocities.  Equation  4 may  be  approximated  by  the  relation 


PH  = plcouo  (5) 

As  stated  in  the  Literature  Survey,  many  researchers  believed 
that  the  shock  pressure  generated  during  the  Impact  of  water  on  a rigid 
surface  could  be  given  by  Equation  5.  Klnslow  (Reference  39)  apparently 
used  this  relation  In  expressing  shock  pressures  for  a water  jet  Impact 
at  633  m/s.  Although  this  relation  may  be  adequate  for  very  low  velocity 
Impact,  It  deviates  markedly  from  Equation  4 for  higher  velocities. 

Figure  2 demonstrates  the  differences  In  pressures  found  from  these  two 
relationships  for  water.  At  Mach  ■ 1.0,  Equation  5 gives  a shock  pressure 
of  only  1/3  that  of  the  actual  value  given  by  Equation  4.  Klnslow's 
results,  for  a velocity  of  633  m/s,  are  seen  to  be  In  error  by  more 
than  40X. 
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Figure  2.  Effect  of  Compressibility  on  the  Hugonlot  Pressure 
for  Water,  (1  MN/m2  ■ 10  Bars) 

b.  Release  Regime 

As  soon  as  the  shock  Is  formed.  It  begins  to  propagate  up  the 
projectile  and  radial  release  waves  propagate  In  towards  the  center  axis 
of  the  projectile.  The  problem  can  no  longer  be  considered  one-dimensional 
In  nature.  For  the  normal  Impact  of  a cylinder,  the  problem  becomes  two- 
dimensional,  axlsymmetrlc. 

Figure  3 shows  the  release  regime  of  impact  for  a cylinder  with 
an  original  L/D  of  2.  Figure  3-b  Illustrates  the  shocked  region  In  the 
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Figure  3.  Shock  end  Release  Waves  In  Fluid  Impact. 

projectile  just  after  Impact,  The  pressure  In  the  shocked  region  Is 
given  by  Equation  4.  It  should  be  noted  that  the  velocities  of  the  shock 
and  release  waves  are  much  greater  than  the  Initial  velocity  of  the 
projectile.  Figure  3-c  shows  the  moment  when  the  release  waves  have 
converged  on  Pt.  8,  the  axis  of  the  cylinder.  Figure  3-d  shows  the 
moment  when  the  release  waves  have  just  caught  the  center  of  the  shock, 
Pt.  C.  The  region  of  fully  shocked  material  no  longer  exists.  The 
curvature  of  the  shock  Is  due  to  the  release  process,  which  has  progres- 
sively weakened  the  shock,  thereby  decreasing  the  shock  velocity. 
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Several  reflections  of  the  release  wave  are  probably  necessary 
to  establish  steady  flow.  The  existence  of  steady  flow,  of  course.  Is 
dependent  on  the  length  (or  L/D)  of  the  cylinder.  For  very  short 
cylinders,  the  Impact  event  will  be  over  before  the  release  regime  has 
ended.  Although  there  Is  no  simple,  analytical  method  for  studying 
the  release  regime,  several  critical  Items  can  be  Identified. 

The  duration  of  the  shock  pressure  at  the  center  of  Impact  can 
be  approximated  by  the  amount  of  time  that  It  takes  the  Initial  release 
wave  to  reach  the  center,  Pt.  B of  Figure  3-c.  The  release  wave  can  be 
considered  to  be  a fan  of  weak  expansion  waves,  which  Implies  that  the 
release  process  can  be  considered  Isentroplc.  The  velocity  of  the 
Initial  release  wave  Is  just  the  speed  of  sound  In  the  shocked  material, 
cr,  which  is  given  by  the  expression 

c2  = (dP/dp  )p  (6) 

r H 

which  is  the  slope  of  the  isentroplc  pressure-density  curve  at  the 
Hugoniot  state.  Thus,  the  expression  for  the  time  necessary  for  the 
release  wave  to  initially  reach  the  center  Is 

tB  * a/cr  (7) 

where  a Is  the  Initial  radius  of  the  cylinder.  Figure  4 shows  the  speed 
of  sound  in  the  shocked  region,  cf,  plotted  versus  Impact  velocity  for 
water.  Using  this  Information,  the  duration  of  the  Hugoniot  pressure, 
tg,  can  be  obtained  from  Equation  7.  Figure  5 shows  the  relationship 
between  tg  and  the  Impact  velocity,  uQ,  for  cylindrical  projectiles  of 
water  of  various  radii.  (The  equation  of  state  of  water  is  given  in  a 
later  section). 

Another  important  item  Is  the  time  that  it  takes  the  front  of  the 
release  wave  to  just  capture  all  of  the  shock  wave.  As  stated  previously, 
immediately  upon  impact  the  release  wave  begins  to  interact  with  the 
shock  at  the  edge  of  the  projectile.  Since  the  wave  speed  In  the 
fully  shocked  medium  is  always  greater  than  the  shock  speed,  the 
release  wave  will  Interact  with  progressively  more  of  the  shock  as  the 
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Figure  4.  Comparison  of  the  Shock  Velocity,  u , and  the  Sound 
Speed  In  the  Shocked  Region,  cr*  5 


Floure  5,  Variation  of  the  Duration  of  the  Hugonlot  Pressure  with 
Impact  Velocity  for  a Cylindrical  Projectile  of  Water. 
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Impact  continues.  Figure  3-d  shows  the  condition  In  which  the  release 
waves  have  just  converged  on  Pt.  C.  After  this  time,  the  pressure  In 
the  region  behind  the  shock  will  rapidly  decay.  The  time  of  this 
occurre-ice,  t£,  may  be  derived  from  geometric  considerations.  In  order 
for  the  release  wave  starting  at  Pt.  A to  Interact  with  the  shock  at 
Pt.  C,  It  must  travel  a radial  distance,  a,  and  an  axial  distance  equal 
to  the  axial  distance  traveled  by  the  shock.  At  the  time  of  Intersection, 
the  shock  has  propagated  a distance 

x.  = (u-u)t  (0) 

S S 0 c 

where  (u$-uo)  Is  the  velocity  of  the  shock  relative  to  the  target  as 
seen  in  Figure  1-b.  From  this  It  can  be  seen  that  the  release  wave  will 
have  traveled  a distance 

x = <x2  ♦ a2)*  (9) 

r s 

In  the  time,  tc,  given  by  the  expression 

tc  = xr/cr  (10) 

Thus,  It  can  be  seen  that 


t 

c 


a 

<c2-(u  -u  )*)H 
r so 


(11) 


From  Equation  11  an  expression  can  be  derived  for  the  critical 
projectile  length,  Lc>  which  is  the  length  for  which  the  radial  release 
wave  will  just  Intersect  the  shock  on  axis,  Pt.  C,  as  the  shock  reaches 
the  end  of  the  projectile.  Again,  from  geometry,  it  can  be  seen  that 


L 


c 


u t 
s c 


(12) 


Combining  Equations  11  and  12  and  nondlmensionallzlng,  It  is  seen  that 


(L/D) 

c 


u 

s 


2{c2-(u  -u  )2}** 
r so 


(13) 


For  a projectile  with  an  L/D  < (L/D)  , the  shock  will  reflect  off  the 
projectile  rear  surface  before  it  has  all  been  captured  by  the  radial 
release  waves.  The  shock  will  be  reflected  In  the  form  of  a rare- 
faction wave.  This  rarefaction  wave  could  complicate  the  impact  event 
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by  causing  a tensile  failure  In  the  projectile,  decreasing  the  Incasing 
velocity  of  the  end  of  the  projectile,  and  disturbing  the  radial  release 
process. 

However,  for  a projectile  with  an  L/D  > (L/D)c . the  shock  will 
be  substantially  weakened  by  the  release  waves  prior  to  reaching  the  pro- 
jectile rear  surface  and  Its  effects  will  be  reduced  or  effectively 
cancelled.  Figure  4 showed  both  u$  end  cf  as  a function  of  Impact 
velocity  for  water.  Using  these  values,  (L/D).  for  a cylinder  of  water 
can  be  obtained  froai  Equation  13.  Figure  6 shows  a plot  of  (L/D)c  versus 
ImpACt  velocity  for  water. 


IMPACT  VELOCITY  , ^ (m/a) 

Figure  6.  Variation  of  Critical  Length  with  Impact  Velocity 
for  Water. 
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For  a projectile  o*  sufficient  length,  steady  flow  should  be  set  up 
up  after  several  reflections  of  the  radial  release  waves.  It  seems  likely 
that  for  a projectile  with  a length  somewhat  greater  than  Lc,  a condition 
of  steady  flow  should  exist  for  some  time. 

c.  Steady  Flow  Regime 

During  the  release  phase,  the  shock  Is  constantly  being  weakened 
by  the  release  waves,  with  an  accompanying  decrease  In  shock  velocity. 

For  a low  velocity  Impact  (subsonic),  the  shock  wave  will  be  weakened  to 
the  point  that  It  will  disappear.  By  this  time,  for  steady  flow,  stream- 
lines will  be  established.  However,  for  a supersonic  Impact,  the  shock 
wave  will  not  disappear.  Instead,  Its  velocity  will  decrease  until  It 
becomes  a standing  shock.  Behind  this  standing  shock,  the  flow  will  be 
subsonic  and  will  follow  steady  state  streamlines.  Because  of  the 
standing  shock,  the  velocity  and  pressure  fields  set  up  In  the  fluid  will 
be  different  for  the  two  cases. 

(1)  Subsonic  Flow 

If  it  Is  assumed  that  the  shearing  forces  and  body  forces 
are  negligible  throughout  the  flow  field  and  that  all  stream  properties 
vary  continuously  in  all  directions,  then  along  each  streamline, 

Bernoulli's  Equation  can  be  written 

f " +/udu  = K (14) 

J P J 


where  K Is  constant  along  the  streamline  and  may  vary  from  one  stream- 
line to  another.  For  the  case  of  a cylinder  impacting  a rigid  plate, 
the  flow  field  is  essentially  uniform  at  some  distance  away  from  the 
impact  surface,  so  that  the  K must  be  the  same  for  each  of  the  stream- 
lines in  this  region.  This  implies  that  K has  the  same  value  throughout 
the  entire  flow  field.  For  this  case,  the  pressure  P,  at  any  point  In 
the  flow  field  can  be  related  to  the  velocity,  u,  at  that  point  by  the 
relation 


(15) 
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where  PQ  and  uQ  are  the  pressure  and  velocity  of  the  uniform  flow  field 
some  distance  away  from  the  Impact  surface  and  are  approximated  by  the 
atmospheric  pressure  and  the  Initial  Impact  velocity. 

In  order  to  obtain  the  pressure  at  any  point  along  the  impact 
surface  from  Equation  15,  the  velocity  at  that  point  must  be  known  and 
the  equation  of  state  of  the  material,  p = p(P),  must  be  known.  The 
equations  of  state  for  soft  body  materials  will  be  discussed  in  a later 
section.  The  expression  for  the  velocity  at  a point  is  not  found  so 
easily.  For  a two-dimensional  fluid  jet,  the  velocity  field  in  the 
fluid  was  obtained  by  Milne-Tbompson  (Reference  55)  from  potential  flow 
theory.  However,  no  closed  form  solution  of  potential  flow  theory  has 
been  developed  for  a circular  jet.  The  normal  approach  is  to  assume  an 
expression  for  the  velocity  field  based  on  empirical  data. 

In  this  approach,  the  pressure  at  the  center  and  the  pressure  at 
the  edge  are  found  using  Equation  15.  Then  a general  expression  for  the 
pressure  distribution  is  assumed  based  on  empirical  data.  This  expression 
is  forced  to  satisfy  the  pressure  boundary  conditions  and  conservation 
of  momentum. 


For  the  normal  impact  of  a uniform,  cylindrical  projectile  on  a 
rigid  plate,  axial  symmetry  dictates  that  the  point  at  the  center  of 
the  plate  be  the  stagnation  point.  The  pressure  at  the  center  of  the 
plate  is  the  stagnation  pressure,  Pg  (gage  pressure),  and  the  velocity 
at  the  center  is  zero.  Therefore,  at  the  center  of  the  plate,  Equation 
15  takes  the  form 


/*(P  8 + VdP 
J p P 


(16) 


If  the  fluid  is  assumed  incompressible.  Equation  16  gives 

po  = h pu  2 


(17) 


For  most  materials,  the  density  tends  to  increase  the  applied  pressure,  so 
that  Equation  16  implies 


p 

s 


* H pu  2 
o 


(18) 
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The  second  pressure  boundary  condition  Is  that  the  pressure  must  go  to 
zero  at  some  radial  distance  from  the  center. 


Simple  momentum  considerations  require  that  during  steady  flow, 
the  Impulse  Imparted  to  the  target  by  a unit  of  fluid  must  be  equal  to 
the  axial  momentum  lost  during  Impact.  This  may  be  written 


(M  du) 


(19) 


where  the  force,  F,  Is  assumed  constant  during  steady  flow  and  the 
duration  of  Impact  is  represented  by  tQ.  For  a unit  of  fluid  with 
Initial  values  of  muss  M,  density  p,  length  L,  and  cross-sectional 
area  A,  this  expression  becomes 

F tD  = M(u-uq)  (29) 

For  most  soft  body  materials,  the  rebound  velocity  after  Impact,  u,  Is 
so  small  that  it  may  be  Ignored.  Also,  during  steady  flow,  the  duration 
of  impact  for  a unit  of  fluid  of  length,  L,  Is  simply  the  time  required 
for  the  fluid  to  flow  through  Its  length.  That  Is, 

tD  = L/uo  <21> 

Thus,  the  force  generated  in  the  steady  flow  regime  is  seen  to  take  the 
form 

F - pAu2  (22) 


*.  Since  the  force  Is  simply  the  Integral  of  the  pressure  over  the  Impact 
surface,  this  expression  may  be  rewritten 

^co 

2tt  / Prdr  = pAu2  (23) 

Jo  ° 

From  the  above  arguments.  It  is  evident  that  any  assumed  expression  for 
pressure  must  predict  the  pressure  to  be  maximum  at  the  center,  given 
by  the  stagnation  pressure  Ps,  decreasing  with  radius  to  a value  of 
zero  at  some  point,  and  that  the  expression  for  pressure  must  satisfy 
Equation  23. 


As  stated  In  the  Literature  Survey,  both  Banks  and  Chandrasekhara 
(Reference  56)  and  Leach  and  Walker  (Reference  58)  developed  expressions 
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for  the  pressure  distribution  due  to  the  normal  Impact  of  a water  jet. 

The  expression  of  Banks  and  Chandrasekhara  took  the  form 

p = **  pu0?-  exp  (24) 

and  that  of  Leach  and  Walker  took  the  form 

t- 

P = >4  PU  * {1-3<J^-)2  -f  (~-)3>  (25) 

where  r Is  the  radial  distance  from  the  center,  a Is  the  Initial  radius 
of  the  jet,  and  ^ and  are  constants  used  to  make  Equations  24  and  25 
conform  to  Equation  23.  For  this  case,  ^ - 0.5  and  c2  “ 2.58. 

Figure  7 shows  a plot  of  Equations  24  ard  25. 

The  main  drawback  to  these  two  expressions  Is  that  they  define  a 

2 

maximum  pressure  at  the  center,  P$,  of  1/2  puQ  . As  was  shown  earlier, 
this  Is  only  true  exactly  for  an  Incompressible  fluid.  For  a compres 
sible  material,  the  pressure  will  be  greater.  The  effect  of  compres- 
sibility on  the  stagnation  pressure  can  be  demonstrated  by 


Figure  7.  Radial  Distribution  of  Steady  Flow  Pressure. 
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considering  the  Impact  of  a cylinder  of  water.  Using  the  equation  of 

state  for  water  which  Is  given  In  a later  section,  Bernoulli's  compressible 

equation  In  the  form  of  Equation  15  can  be  solved  for  P^.  Figure  8~a 

shows  the  results.  Although  the  stagnation  pressure  Is  not  dramatically 

2 

different  from  1/2  puQ  for  low  velocities,  the  difference  Is  evident 
at  higher  velocities.  This  pressure  difference  would  be  even  greater 
for  very  compressible  materials. 


For  soft  body  impacts.  Equations  24  and  25  may  be  generalized  to 
allow  for  the  Increased  stagnation  pressure  due  to  compressibility.  Thus, 
they  may  be  written 


and 


P * Ps  txp 


P . ps  ♦ 


(26) 


(27) 


respectively,  where  the  constants  c-j  and  ;2»  derived  In  Appendix  A,  are 
given  by  Equations  A-7  and  A-9 


C2  = (3.33  Pu^/Pg)*5 

The  effect  of  compressibility  on  the  pressure  distribution  should 

be  evident.  Figure  8-a  shows  that  for  an  impact  velocity  of  uA  a 1400  m/s, 

2 0 
P$  “ 1.14  (l/2ouQ  ).  Figure  8-b  shows  a plot  of  Equation  27  using  this 

value  of  P . 

(2)  supersonic  Flow 

For  the  supersonic  Impingement  of  a fluid  jet  normally 
onto  a rigid  plate,  a standing  shock  Is  set  up  In  the  flow.  The  change 
in  particle  velocity  across  the  shock  may  be  small.  However,  the  flow 
behind  the  shock  will  be  subsonic  because  of  the  greatly  Increased 
local  sound  speed. 
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Figure  8.  Compressibility  Effects  In  Steady  Flow,  a)  Stagflation 
Pressure,  b)  Pressure  Distribution. 
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In  the  previous  section  on  subsonic  flow,  Bernoulli's  Equation 
was  used  to  relate  the  pressures  along  the  Impact  surface  to  the  pressure, 
PQ,  density  pQ,  and  velocity,  uQ,  In  the  uniform  flow  region.  However, 
for  supersonic  flow  this  approach  cannot  be  used  since  the  standing  shock 
Is  a discontinuity  surface  and  Bernoulli's  Equation  only  holds  In  a 
continuous  flow  field.  Instead,  the  change. In  stream  properties  across 
the  shock  must  be  evaluated  using  shock  relations.  Bernoulli's  Equation 
can  then  be  used  throughout  the  flow  field  behind  the  shock.  Figure  9 
Illustrates  this  case.  Region  1 In  front  of  the  shock  Is  a uniform  flow 
field  and  the  stream  properties  In  this  region  are  the  Initial  properties, 
namely  PQ,  pQ,  and  uQ.  In  Region  2 just  behind  the  shock,  the  properties 
may  not  be  constant  over  the  cross  section  since  the  shock  may  not  be  a 
perfectly  plane  shock.  However,  at  the  center  axis,  the  shock  can  be 
assumed  plane  and  one-dimensional  shock  equations  can  be  used.  The 
properties  along  the  axis  just  behind  the  shock  Will  be  designated 
Pg,  P2»  U2« 


Figure  9.  Standing  Shock  In  Steady  Flow. 
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Since  no  uniform  flow  field  exists  behind  the  shock.  Bernoulli's 
Equation  will  take  the  form 


(28) 


and  will  hold  only  along  the  streamline  which  passes  through  the  center 
of  the  projectile  and  the  stagnation  point  In  the  center  of  the  Impact 
region. 


The  case  of  a standing  shock  was  shown  In  Figure  1-c.  For  this 


case.  Equations  4 and  1 take  the  form 


p2  * poVp 


(29) 


where 


poUs  * p2(Vup) 


u = u 
s o 


(30) 

(31) 


u = u -u_ 
p o 2 


(32) 


One  more  equation  Is  needed  in  order  to  solve  this  set  of  equations. 


Such  an  equation,  relating  the  shock  velocity,  us,  with  the  particle 


velocity,  up,  will  be  given  for  various  soft  body  materials  In  a later 


section.  Thus,  for  a given  set  of  Initial  conditions  and  a given  material. 
Equations  29,  30,  31,  and  32  can  be  solved  for  Pg,  Pg,  and  Ug. 


Using  these  values,  the  pressure  at  the  stagnation  point  on  the 
surface  of  the  plate  can  be  found  using  Bernoulli's  Equation  of  the  form 


r*^.  v 
j Pj  ” — 


(33) 


In  this  relationship,  the  equation  of  state  of  material  must  be  the 
pressure-density  relationship  for  an  Isentroplc  compression  passing  through 
the  points  P?,  p?. 
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As  In  the  subsonic  flow  case,  once  P$  Is  determined,  Equations  26 
and  A-7  or  27  and  A-9  can  be  used  to  develop  an  expression  for  the 
pressure  distribution  along  the  surface  of  the  plate. 

d.  Termination  of  Impact 

As  stated  previously,  the  particles  of  fluid  flow  along  stream- 
lines which  are  set  in  space.  These  streamlines  are  curved  near  the 
plate  surface  so  that  the  Impacting  material  is  "turned"  during  Impact. 

For  subsonic  flow,  as  the  fluid  nears  the  target  surface  there  is  a 
gradual  decrease  in  velocity  with  a corresponding  Increase  in  local  pres- 
sure. Thus,  during  steady  flow  there  Is  a pressure  field  set  up  In  the 
fluid  In  which  the  local  pressure  is  maximum  at  the  stagnation  point  and 
decreases  to  the  atmospheric  pressure,  P , at  a substantial  distance 
from  this  point.  As  the  end  of  the  projectile  nears  the  impact  surface, 
it  enters  this  field  of  increasing  local  pressure  and  disrupts  the  field 
due  to  the  immediate  drop  in  pressure  behind  the  projectile  rear  surface. 
For  supersonic  impact,  this  happens  when  the  projectile  end  passes 
through  the  standing  shock.  In  either  case,  release  waves  emanate  from 
the  re'r  surface  and  propagate  to  the  impact  surface,  thereby  causing  a 
slin^t  decrease  in  the  impact  velocity  of  the  rear  surface  as  well  as  a 
d.  se  in  t pressure  along  the  impact  surface  of  the  target.  This 
p.  ss  continues  until  the  end  of  the  projectile  reaches  the  surface  of 
the  plate  and  the  impact  event  is  ended. 

As  stated  earlier  in  Equation  21,  the  total  duration  of  the 
impact  can  be  approximated  by  the  time  needed  for  the  projectile  to 
"flow  through"  its  length,  or  tQ  3 l/uQ.  Deviation  from  this  could  be 
caused  by  the  decrease  in  projectile  velocity  due  to  release  waves  from 
the  back  of  'ojectlle. 

2.  SPECIAL  CONSIDERATIONS 

In  the  last  section  a theory  was  presented  concerning  the  normal 
impact  of  a right-circu  cylindrical  projectile  of  fluid  against  a 
rigid  target.  However,  there  are  several  other  cases  of  special  interest. 
Included  among  these  is  the  oblique  impact  of  a right-circular  cylinder, 
the  normal  impact  of  a yawed  projectile,  and  the  normal  impact  of  a 
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projectile  with  a curved  leading  edge.  Next,  the  porosity  of  the 
projectile  material  Is  seen  to  have  a great  effect  on  the  impact  pres- 
sure. Finally,  the  Impact  of  a fluid  projectile  against  a nonrlgld 
target  Is  considered. 

a.  Projectile  Yaw 

Figure  10-a  presents  the  normal  Impact  of  a right-circular 
cylinder  In  which  the  axis  of  the  cylinder  Is  yawed  by  an  angle  $ to  the 
velocity  vector.  The  amplitude  of  the  shock  Hugonlot  pressure  will  be  the 
same  as  that  for  a projectile  Impacting  without  yaw.  However,  the  dura- 
tion of  this  pressure  at  various  points  along  the  target  surface  will 
differ  If  the  yaw  angle,  <J>,  is  greater  than  the  critical  angle,  <{>cr. 


Figure  10.  The  Impact  of  a Yawed  Projectile. 


In  Figure  10-b,  the  projectile  Is  seen  during  Initial  impact. 
If  the  shock  originating  at  point  a arrives  at  point  b prior  to  point  b 
reaching  the  impact  surface,  then  a release  wave  will  develop  and  propa 
gate  back  Into  the  shocked  region.  The  critical  angle,  <Pcr . is  the 
angle  of  yaw  for  which  the  shock  will  just  reach  point  b at  the  Instant 
point  b reaches  the  surface.  From  the  figure,  it  can  be  seen  that  for 
<j>  * 4^,  the  ratio  of  the  distance  traveled  by  point  b,  be,  to  the 


distance  traveled  by  the  shock,  ab,  must  be  equal  to  the  ratio  of  the 
velocity  of  point  b,  u.  to  the  shock  velocity,  u . This  may  be  written 

O « 


for  A = A . 
T Tor 


From  geometry,  It  Is  evident  that 


sin  $ 

cr 


(34) 


(35) 


Thus,  the  relationship  for  the  critical  angle  becomes 


(36) 


The  critical  angle  changes  with  velocity  for  most  materials  since  the 
shock  velocity  Is  normally  a function  of  the  Impact  velocity. 

For  an  impact  In  which  the  yaw  is  less  than  the  critical  angle, 
the  amplitude  and  duration  of  the  shock  pressure  at  the  various  points 
over  the  surface  will  be  essentially  the  same  as  that  for  the  impact  of 
a projectile  without  yaw.  However,  for  an  impact  In  which  the  yaw  Is 
greater  than  the  critical  ancle,  the  full  shock  Hugoniot  pressure  will 
exist  only  at  the  initial  point  of  impact.  At  each  of  the  other  points 
on  the  impacting  surface  of  the  projectile,  the  material  will  have  been 
shocked  and  partially  released  before  Impacting  the  surface,  thereby 
reducing  the  shock  pressure  generated  upon  impact.  Also,  the  duration 
of  the  shock  pressures  at  the  center  will  be  decreased  due  to  the 
immediate  formation  of  release  waves. 

b.  Projectile  Leading  Edge  Curvature 

For  the  Impact  of  a curved  surface  on  a flat  target,  the  angle, 
4>,  which  the  Impacting  surface  makes  with  the  target  surface  Increases 
as  the  impact  progresses.  This  condition  is  demonstrated  ir>  Figure  11. 
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Figure  11.  The  Impact  of  a Curved  Surface  on  a Flat  Target. 

Once  the  angle  <p  exceeds  $cr  as  def1ned  by  Equation  36,  a release  wave 
will  form  and  propagate  Into  the  shocked  region.  This  Implies  that  the 
duration  of  the  Hugonlot  shock  at  the  center  will  be  less  for  the  Impact 
of  a sphere  than  for  a cylinder  <w1thout  yaw)  of  like  diameters s 

c.  Oblique  Impact 

Figure  12-a  demonstrates  the  oblique  Impact  of  a right-circular 
cylinder  against  a rigid  plate  In  which  the  velocity  vector  Is  directed 
along  the  axis  of  the  cylinder  and  Intersects  the  target  at  an  angle  a 
relative  to  the  target  surface. 

(1)  Shock  Regime 

For  the  oblique  Impact  of  a projectile  on  a rigid  plate, 
a coordinate  transformation  aids  in  the  understanding  of  the  shock 
process.  From  Figure  12-a  it  Is  evident  that  the  component  of  projectile 
velocity  normal  to  the  plate  Is  uQ  sin  a and  the  component  tangential 
to  the  plate  Is  uQ  cos  a.  Figure  12-b  shows  the  impact  in  the  transformed 


Figure  12.  Oblique  Impact  of  a Cylinder  on  a Rigid  Plate. 
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coordinate  system  In  which  the  projectile  is  seen  to  impact  normally  onto 
the  moving  plate.  Assuming  the  frictional  forces  along  the  surface  are 
, negligible,  this  case  Is  identical  to  the  normal  Impact  of  a yawed 
projectile  with  an  Initial' velocity  of  uQ  sin  a.  Correspondingly,  the 
Hugonlot  shock  pressure  will:  be  Identical  to  that  produced  by  a normal 
Impact  of  a projectile  with  an  Initial  velocity  of  uQ  sin  a.  Figure  13 
demonstrates  the  decrease  In  shock  pressure  with  Impact  angle  of  obliquity 
for  the  Impact  of  a cylinder  of  water. 

/ 

(2)  Release  Regime  ! 

The  release  process  Is  similar  to  that  of  the  yawed 
projectile.  For  this  case,  the  angle  that  the  impacting  projectile 
surface  makes  with  the  target  surface  Is  seen  to  be  governed  by  the 
relation  $ * 90°-a.  For  the  Impact  In  which  <{>  > <frcf,  the  full  Hugonlot 
pressure  will  occur  only  at  the  initial  point  of  Impact. 

(3)  Steady  Flow 

The  steady  flow  phase  of  Impact  for  a fluid  cylinder  is 
shown  in  Figure  14.  From  momentum  considerations  it  can  be  seen  that 
the  majority  of  fluid  will  flow  "downstream".  This  causes  the  stagnation 
point  to  shift  "upstream"  of  the  center  of  the  plate.  As  long  as  a 
stagnation  point  exists,  the  pressure  at  that  point  will  be  given  by 
Equation  16.  Thus,  the  maximum  pressure  generated  during  steady  flow 
will  be  Independent  of  the  angle  of  Impact.  However,  the  distribution 
of  pressure  along  the  surface  will  be  greatly  dependent  on  the  impact 
angle. 


Although  no  expression  for  pressure  distribution  in  an 
oblique  impact  of  a circular  jet  has  been  developed,  Taylor  (Reference  59) 
developed  a distribution  for  a two-dimensional  (plane)  jet  of  water.  This 
distribution  Is  shown  in  Figure  15  for  an  impact  angle,  a,  of  30°. 

In  this  figure,  the  profile  of  the  projectile  Is  superimposed  on  the 
pressure  distribution  tc  demonstrate  the  distance  that  the  stagnation 
point,  s,  has  shifted  upstream  from  the  center  of  impact. 
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Figure  14.  The  Steady  Flow  Phase  of  an  Oblique  Impact. 
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Figure  15.  Steady  Flow  Pressure  Distribution  for  a 30°  Impact  of 
a Plane  Jet  of  Water  (Taylor,  Reference  59). 

From  simple  momentum  considerations.  Equation  22  stated  that 
the  total  force  exerted  on  the  target  at  any  Instant  of  time  for  a normal 

p 

impact  could  be  given  by  the  relation  F * pAuQ . However,  the  momentum 
transferred  to  the  target  during  an  oblique  Impact  is  only  the  normal 
component  of  the  total  momentum.  Thus,  for  an  oblique  impact.  Equation  22 
can  be  rewritten 

F = pAu*  sin  a (37) 

For  normal  impact,  a = 90°  and  Equation  37  reduces  to  Equation  22.  Thus, 
the  pressure  distribution  for  an  oblique  impact  is  seen  to  be  different 
from  that  for  a normal  impact,  due  both  to  the  decrease  in  total  force 
and  the  loss  of  axial -symmetry.  , 

i * 

t 

d.  Material  Porosity 

In  the  previous  section  on  Hydrodynamic  Theory,  the  Hugoniot 
shock  pressure  was  presented  as  a function  of  Impact  velocity  for  water 
in  Figure  2.  Later,  in  Figure  7,  the  distribution  of  steady  flow  pres- 
sure was  given  for  the  impact  of  a jet.  of  water.  In  both  of  these  cases, 
the  water  was  assumed  to  have  no  porosity.  As  will  be  seen  In  the  next 
section  on  Equations  of  State  of  soft  body  materials,  the  presence  of 


I 
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porosity  has  a great  effect  on  both  the  shock  velocity  and  compressibility 
of  a material  during  Impact.  Thus,  the  equation  of  state  for  a porous 
material  differs  from  that  of  a nonporous  material. 

The  decrease  In  shock  velocity  due  to  prooslty  results  In  a 
marked  decrease  In  the  shock  pressures.  Figure  16  presents  the  variation 
In  Hugonlot  pressure  with  porosity,  z,  for  the  normal  impact  of  a cylinder 
of  water.  (As  defined  In  the  next  section,  z Is  simply  the  volume 
fraction  of  air  In  the  water).  Likewise,  the  decrease  In  density  due 
to  the  addition  of  porosity  causes  a decrease  In  the  stagnation  pressure 
during  steady  flow.  However,  the  Increase  In  compressibility  effects 
tend  to  counteract  this,  so  that  the  resulting  decrease  Is  relatively 
small.  Figure  17  presents  the  variation  In  the  steady  flow  stagnation 
pressure  for  water  with  porosity. 

In  Figure  8,  the  alteration  of  the  steady  flow  pressure 

distribution  due  to  compressibility  effects  was  presented.  In  that 

case,  the  Increased  material  compressibility  resulted  from  Increased 

Impact  velocity.  Compressibility  effects  are  even  more  pronounced 

for  materials  which  have  porosity.  Figure  18  presents  the  pressure 

distribution  for  both  porous  and  nonporous  water  using  Equation  27. 

In  this  figure,  the  pressure  is  nondlmotislonallzed  with  respect  to 

2 

the  Incompressible  flow  pressure  l/2pzu0  , where  p„,  defined  In  the  next 
section  by  Equation  58,  Is  the  average  density  of  the  porous  water  for 
this  case.  Also,  the  radial  distance  from  center  Is  nondimenslonalized 
relative  to  the  Initial  projectile  radius,  a. 

e.  Non-Rigid  Target 

Up  to  this  point  the  theory  has  been  concerned  with  impact  on  a 
rigid  plate.  However,  all  real  target  materials  are  non-Hg1d  to  son* 
extent.  Thus,  the  response  of  the  target  during  impact  cannot  be 
overlooked.  If  the  target  material  can  be  assumed  to  behave  elastically, 
the  Initial  target  response  is  straightforward.  At  very  early  times 
during  the  Impact,  only  the  local  area  of  the  target  directly  under 
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Figure  16.  Relationships  Between  the  Hugoniot  Shock  Pressure  and 
Particle  Velocity  for  Water  with  Various  Porosities. 

(1  MN/m2  « 10  Bars) 


Figure  17.  Relationships  Between  the  Stagnation  Pressure  and  the 
Impact  Velocity  In  Steady  Flow  for  Water  with  Various 

"’oroslties.  (1  MN/to2  “ 10  Bars) 


42 


AFML-TR-77-134 


Figure  18.  The  Variation  In  Steady  Flow  Pressure  Distribution 
for  Water  Due  to  the  Addition  of  Porosity. 


the  Impacting  projectile  Is  affected  by  the  Impact.  Figure  19  illustrates 
this  early  shock  regime.  From  this  figure  It  is  evident  that  to  a first 
approximation  the  shock  waves  may  be  considered  Initially  planar,  and 
ons-dlmenslonal  shock  equations  can  be  used  for  the  target  as  well  as 
the  projectile. 
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Figure  19.  impact  on  an  Elastic  Target  During  the  Early 
Shock  Phase. 
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In  Figure  19,  Region  2 Is  the  shocked  region  In  the  projectile 
and  Region  3 Is  the  shocked  region  In  the  target. 

For  the  Impact  on  a non-rigid  target,  Equation  4 Is  seen  to 
take  the  form 

ph  = Pusup  08) 


where  up,  the  particle  velocity,  Is  no  longer  equal  to  the  Initial  Impact 
velocity.  In  Regions  2 and  3,  Equation  38  may  be  written 


P2  = pp  us  Un 

P V 


p,  = p.  uo  u 
3 t st  pt 


(39) 

(40) 


At  the  projectile- target  Interface,  equilibrium  requires  that 


p = p 
2 r3 


U2  * U3 


From  the  definition  of  particle  velocity,  It  can  be  seen  that 


(41) 


u = u -u* 
Fp  o 2 


(42) 


V = U3'U4  = U3 


(43) 


Thus,  the  particle  velocities  can  be  related  by  the  expression 


u = u -u 
p*  '■>  F, 


(44) 


Substituting  Equation  44  Into  Equation  40  and  eouating  the  pressures  F2 

p u u - p .u  (u  -u  ) (4i>) 

p SP  PP  t PP 


and  gives 
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From  this,  it  can  be  seen  that 


ptV 


p U + p . u 

p Sp  t Bt 


(46) 


Thus,  the  pressures  generated  at  the  projectile-target  interface  for 
the  impact  of  an  elastic  target  material  are  seen  to  be 


p 


p u u 
p s o 
* p 


•VV 


•vy  p*\ 


(47) 


To  demonstrate  the  difference  between  the  pressures  predicted 
by  Equation  4,  which  assumes  a rigid  target,  and  Equation  47,  several 
calculations  were  made  for  the  normal  Impact  of  water  on  various  targets 
at  a velocity  of  200  m/s.  For  the  Impact  on  steel,  the  shock  pressure 
found  using  Equation  47  is  approximately  4%  less  than  that  found  using 
Equation  4.  For  titanium,  the  difference  is  approximately  8%;  for 
aluminum,  approximately  115!;  and  for  polycarbonate,  approximately  355!. 

For  thin  targets,  the  Initial  shock  wave  in  the  target  reflects 
off  the  rear  surface  as  a tensile  wave  of  similar  strength.  This  wave 
propagates  back  to  the  impact  surface,  causing  a decrease  in  the  shock 
pressures  at  the  surface.  Continuing  wave  reflections  cause  an  increase 
in  the  local  particle  velocity  of  the  target  and  a decrease  in  pressure. 
As  the  impact  proceeds,  the  entire  target  structure  develops  a velocity 
along  the  initial  axis  of  impact.  The  net  result  of  this  target 
deformation  is  that  the  relative  velocity  between  the  projectile  and 
target  decreases  with  a resulting  decrease  in  the  flow  pressures. 


3.  EQUATIONS  OF  STATE 

In  order  to  determine  the  pressures  generated  during  impact,  the 
material  properties  of  the  projectile  and  target  must  be  known.  These 
are  applied  to  the  problem  through  the  equation  of  state  of  the  materials. 
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a.  Shock  Compression 

The  expression  for  the  Initial  shock  pressure  is  given  by 
Equation  4,  PH  = pQusu0.  In  this  expression,  pQ,  the  initial  density, 
and  uQ,  the  Initial  velocity,  are  known  for  any  given  impact.  Thus,  in 
order  to  solve  for  the  shock  pressure,  a second  expression  is  needed 
whicl  ' js  the  response  of  the  projectile  material  under  shock  loading. 

It  has  been  found  that  for  most  solids  and  fluids  (including 
wc.;er  and  air)  the  relationship  between  the  shock  velocity  and  particle 
velocity  can  be  expressed  as 


u = c+  k u 

so  p 


(48) 


where  k is  a constant  for  the  material  and  cQ  is  the  sound  speed 
(velocity  of  propagation  of  an  infinitesimal  disturbance)  in  the 
material.  This  relationship,  often  called  the  "linear  Hugoniot",  will 
be  usea  to  represent  all  soft  body  materials  in  this  research.  It  can 
be  substituted  into  Equation  4 directly  (by  recalling  that  for  normal 
impact,  of  a rigid  target,  the  particle  velocity,  up,  is  simply  the 
initial  velocity  of  the  projectile).  Thus,  Equations  4 and  48  are 
sufficient  for  determining  initial  shock  pressure. 

However,  it  is  also  desirable  to  obtain  a relationship  between 
the  pressure  and  density  on  both  sides  of  the  shock  by  combining 
Equation  48  with  Equations  1 and  4.  Equation  1,  the  continuity  equation 
across  the  shock,  may  be  written  in  the  form 


u 


-£  = i - 


= q 


(49) 


where  the  subscripts  1 and  2 refer  to  Region  1 in  front  of  the  shock 
and  Region  2 behind  the  shock. 
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Substituting  this  relationship  Into  Equation  48  gives 


ua/cQ  = 1/ (l-kq) 


(50) 


and  Equation  4 may  be  written  In  the  form 


P2  = PlVp 

Thus,  substituting  Equations  49  and  50  Into  Equation  4 gives 


P2  = plCo 


J^(l- 


kq)5 


(51) 


(52) 


Another  approximate  expression  for  the  pressure-density  relation- 
ship across  a one-dimensional  shock  was  given  by  Cogolev,  et  al. 

(Reference  60),  and  takes  the  fornl 


P2 


fp2  S 


(53) 


where  A and  B are  material  constants.  Ruoff  (Reference  61)  has  demon- 
strated that  for  a material  which  exhibits  a linear  Hugonlot,  Equation  48, 
the  constants  may  be  approximated  by  the  expressions 

A = p.jC^/nk-l) 

* * _ - 


where  k is  the  constant  in  Equation  48. 

Thus,  the  pressure  and  density  behind  the  shock  can  be  obtained 
from  either  Equations  52  or  53  if  the  relationship  between  the  shock 
velocity  and  particle  velocity  can  be  approximated  by  Equation  48. 

For  water,  Heymann  (Reference  35)  showed  that  Equation  48  provides  a 
quite  accurate  fit  to  the  data  up  to  Mach  1.2  (approximately  1800  m/s). 
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b.  Isentroplc  Compression 

The  expression  for  the  steady  flow  pressure  is  Bernoulli's  com- 
pressible equation.  Equation  15, 

r ^ r 

I - - + I udu  = 0 

J?  P J u 

o o 


In  order  to  solve  this  equation  ?'ur  a given  impact  velocity,  the 
pressure-density  relationship  v,;  isentroplc  compression  must  be  given 
for  the  projectile  material. 

For  vary  low  pressure  levels,  two  simple  expressions  are  often 
used  to  approximate  the  Isentropic  pressure-density  relationship  for 
solids  and  liquids.  The  first  is  simply 

p = constant  (55) 

Although  never  exact,  the  assumption  of  incompressibility  can  often  be 
used  for  low  velocity  impacts  without  causing  considerable  error. 

This,  of  course,  will  not  be  true  for  porous  materials. 


The  second  relationship  is  obtained  from  the  definition  of  the 
bulk  modulus 


6 = -P 


dP 

dp 


(56) 


For  very  low  pressures,  3 , the  bulk  modulus,  is  often  assumed  constant. 
For  more  accuracy,  a relationship  of  the  form 

3 = Bo  + + $2P2  + (57) 

may  be  assumed. 

However,  for  the  pressures  normally  generated  during  soft  body 

Impacts,  the  assumptions  of  incompressibility  and  constant  bulk  modulus 

are  inadequate.  Instead,  equations  of  state  accurate  to  pressures  of 

o 

several  hundred  MN/m  (several  kilobars)  are  needed.  In  this  pressure 
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range,  the  Isentropic  pressure-density  relationships  for  most  solids 
and  fluids,  Including  water,  can  be  approximated  by  the  Hugonlot  shock 
relationship  (Reference  62)  given  by  Equations  52  or  53. 

c.  Effect  of  Porosity 

It  has  been  stated  that  for  most  soft  body  materials.  Equations 
52  or  53  can  be  used  to  describe  the  pressures  and  densities  resulting 
from  both  Isentropic  compression  and  shock  compression.  Both  of  these 
equations  are  based  on  a linear  Hugonlot  relationship,  Equation  48. 

It  has  been  observed  that  departure  from  linearity  Is  usually  due  to 
porosity  or  phase  change  (Reference  63). 

Although  phase  changes  do  not  occur  In  the  soft  body  Impacts  of 
interest,  many  of  the  soft  body  materials  are  porous.  The  cellular 
structure  of  animal  and  vegetable  matter  normally  contains  a small 
amount  of  porosity.  Rubber,  in  Its  many  uses,  is  often  used  with  large 
amounts  of  porosity.  Thus,  for  these  materials.  Equations  52  and  53  are 
not  representative. 

In  order  to  develop  Isentropic  and  shock  relationships  which 
will  hold  for  porous  materials,  a simple  theory  by  Torvlk  (Reference  64) 
will  be  considered.  Torvik's  theory  predicts  the  pressure-density 
relationship  across  a shock  for  a homogeneous  mixture  based  on  the 
pressure-density  relationships  of  each  constituent.  Although  Torvlk 
used  this  theory  to  develop  a pressure-density  relationship  under  shock 
compression,  there  Is  no  reason  why  It  should  not  be  used  also  to 
develop  the  isentropic  relationship.  Thus,  in  this  report,  Torvik's 
theory  has  been  used  to  develop  both  shock  and  Isentropic  relationships 
for  mixtures  of  which  one  component  <s  air.  The  derivation  of  both  the 
shock  and  isentropic  relationships  using  this  theory  are  presented  in 
Appendix  B.  The  main  assumptions  made  are  that  the  porous  material  is 
macroscopically  homogeneous  and  isotropic  and  that  the  density  of  each 
constituent  of  the  mixture  is  the  same  as  that  for  a homogeneous  sample 
of  that  constituent  of  the  same  pressure. 
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In  the  context  of  this  report,  a porous  material  Is  defined 
simply  as  a mixture  In  which  one  of  the  constituents  Is  air.  The 
amount  of  porosity  of  the  material,  z,  will  be  defined  as  the  volume 
fraction  of  air  In  the  mixture.  If  the  full  density  of  the  soft  body 
material  (density  with  no  porosity  present)  is  represented  by 
and  the  average  density  of  the  soft  body  material  with  porosity  Is 
represented  by  pz>  then  from  Equation  B-l  It  can  be  seen  that  the 
average  density  of  the  material  can  be  defined  by  the  expression 

p2  = ? Pair  + (l-z)  Pf  (58) 

where  the  volume  fraction  of  air,  fvn,  has  been  replaced  by  z and 
the  volume  fraction  of  the  soft  body  material  by  the  term  (1-z). 

For  the  materials  of  Interest  p^»  Pa^rf  so  that  Equation  58  can  be 
approximated  by  the  relationship 

p = (l-/.)  Pf  ^59) 

Z r 

For  a porous  soft  body  material,  the  isentroplc  pressure-density 
relationship  can  be  obtained  from  Equation  B-9.  Substituting  the  terms 
just  defined,  this  equation  becomes 


where  the  subscript  1 refers  to  the  Initial  (unstressed)  state  and  the 
subscript  2 refers  to  the  final  (stressed)  state,  A *nd  B are  empirical 
constants  given  by  Equation  54,  and  a is  the  ratio  of  specific  heats 
of  air.  Similarly,  the  shock  pressure-density  relationship  can  be 
obtained  from  Equation  B-l 8 and  takes  the  form 


-=>  (r  + 1) 


+ z (1-q) 


..Ljud 
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where  q Is  defined  In  Equation  B-17.  The  difference  In  the  Isentroplc 
and  shock  relations  Is  due  to  the  fact  that  the  pressure-density 
relationships  for  air  are  different  for  these  two  cases,  whereas  these 
relationships  have  been  assumed  the  same  for  the  soft  body  material. 

The  effect  of  porosity  can  best  be  demonstrated  by  applying 
Equations  60  and  61  to  the  case  of  water  with  varying  amounts  of  porosity 
added.  The  properties  of  water  are  given  in  Appendix  C.  Figure  20  shows 
two  plots  of  Equation  61  for  this  case  In  which  porosity  Is  varied  from 
0 to  40*.  Figure  20-a  shows  the  variation  In  specific  volume  (density"1) 

O 

with  pressure  up  to  500  MN/m  (5  kbar)  and  Figure  20-b  looks  at  the  lower 

2 

portion  of  these  curves  for  pressures  up  to  20  MN/m  . The'  relationship 
between  the  particle  velocity  and  the  shock  velocity  for  porous  materials 
Is  obtained  by  simultaneous  solution  of  Equations  1.  3.  and  61.  This 
relationship  Is  shown  in  Figure  21  for  water  with  polity  added. 

From  the  solution  of  these  three  equations,  the  relationship  between  the 
shock  pressure  and  particle  velocity  Is  also  obtained.  This  relationship 
was  used  in  an  earlier  section  to  determine  the  effect  of  porosity  on  the 
shock  pressures  generated  during  Impact.  Figure  16  shows  the  Hugoniot 
shock  pressure  as  a function  cf  impact  velocity  and  porosity  for  the 
impact  of  water  onto  a rigid  plate. 

The  isentropic  pressure-density  relationship.  Equation  60,  for 
water  with  porosity  is  shown  in  Figure  22  (for  water,  volume  ■ (density-1). 
In  isentropic  compression,  the  volume  of  air  in  the  mixture  is  seen  to 
vanish  under  substantial  pressure.  Simultaneous  solution  of  Equations  15 
and  60  result  in  the  relationship  between  steady  flow  pressure  and  impact 
velocity.  This  relationship  was  used  In  a previous  section  to  derive 
the  variation  in  the  steady  flow  pressure  at  the  stagnation  point  to 
the  amount  of  porosity  in  a projectile  of  water.  This  information  was 
presented  in  Figure  17. 
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Figure  20.  Hugoniot  Pressure- Volume  Relationships  for  Water  with 
Various  Porosities.  (1  MN/m2  = 10  Bars) 


SHOCK 


Figure  21.  Relationships  between  the  Shock  Velocity  and  the  1 

Particle  Velocity  for  Water  with  Various  Porosities. 
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Figure  22.  Isentroplc  Pressure- Volume  Relationships  for  water 
with  Various  Porosities.  (MN/m2  * 10  Bars) 
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4.  LIMITATIONS  OF  THE  THEORY 

The  basic  theory  assumes  that  a soft  body  can  be  treated  as  a 
homogeneous  fluid  during  impact.  The  main  case  described  was  the 
impact  of  a right-circular  cylinder  against  a rigid  plate,  with  other 

‘ ’ . • r . 

cases  being  studied  briefly.  The  extent  to  which  the  theory  Is  limited 
in  application  will  now  be  considered. 

Application  of  this  theory  to  nonhomogeneous  materials  may  result 
in  several  errors.  Variations  of  density  in  the  material  will  result 
in  variations  of  the  flow  pressure.  Secondly,  mixture  theory,  which 
is  the  basis  of  the  equations  of  state  for  porous  materials,  requires 
that  the  material  be  macroscopically  homogeneous.  If  this  is  not  so, 
the  propagation  of  waves  in  the  material  will  be  affected,  wi*h  an 
accompanying  deviation  in  the  shock  velocity  and  shock  pressure. 

As  stated  earlier,  a material  may  be  treated  as  a fluid  du-ing 
impact  if  the  pressures  generated  during  impact  are  much  greater  than 
the  strength  of  the  material.  For  a projectile  with  a length-to-diameter 
ratio  (L/D)  of  approximately  one  or  less,  this  probably  requires  that  the 
Hugoniot  shock  pressure  be  an  order  of  magnitude  greater  than  the  tensile 
strength  or  shear  strength  of  the  material.  For  a projectile  with  an 
L/D  greater  than  one,  there  is  an  added  requirement  that  the  stresses 
generated  in  the  projectile  during  steady  flow  must  also  exceed  the 
material  strengths.  However,  as  will  be  discussed  later  in  the  Results 
section,  it  appears  that  the  theory  may  accurately  predict  pressure  for 
impacts  in  which  there  are  large  distortions  during  impact,  even  though 
the  material  strength  is  not  exceeded. 

The  theory,  originally  developed  for  a cylindrical  projectile,  is 
not  limited  to  this  geometry.  However,  it  is  probably  the  most  convenient 
to  study,  noth  theoretically  and  experimentally.  For  projectiles  which 
are  spherical  in  shape,  there  is  not  sufficient  length  for  stead.v-flow  to 
be  established.  However,  the  initial  shock  and  release  processes  as 
explained  in  the  theory  should  still  exist.  In  order  for  steady-f’ow  to 
be  established,  the  projectile  must  have  an  L/D  greater  than  one. 
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Finally,  the  theory  requires  that  the  target  remain  flat,  although 
limited  deflections  of  the  target  are  allowed.  Large  distortions  of  the 
surface,  such  as  "pocketing"  of  a plate,  are  not  allowed  since  the  paths 
of  the  streamlines  will  be  altered,  thereby  Influencing  the  flow  pres- 
sure. Also,  the  aerial  dimensions  of  the  target  surface  should  be 
several  times  the  diameter  of  the  projectile,  so  that  full  turning  of 
the  projectile  results. 
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SECTION  IV 

EXPERIMENTAL  TECHNIQUE 

The  purpose  of  the  experimental  program  was  to  measure  the  forces 
and  pressures  generated  during  soft  body  Impact  on  a rigid  surface.  Two 
basic  tests  were  used  In  this  program.  A Hopklnson  bar  with  strain  gages 
was  used  to  measure  Impulse  Imparted  to  the  target  during  the  Impact  and 
a flat  plate  with  pressure  transducers  was  used  to  measure  the  temporal 
distribution  of  pressure  at  various  points  on  the  plate  surface. 

The  Hopklnson  bar  gave  a reliable  measure  of  the  Impulse  (or  momentum) 
Imparted  to  the  Impact  surface  as  well  as  rhe  relative  distribution  of 
imoulse  during  the  Impact  event.  The  flat  plate  gave  a reliable  measure 
of  pressure  Imparted  to  the  target  during  the  Impact  although  there 
were  an  insufficient  number  of  pressure  transducers  to  give  the  distri- 
bution of  pressures  over  the  entire  surface  for  any  given  test.  Both  of 
these  tests  will  be  presented  in  greater  detail  along  with  the  gun 
system  used  for  launching  the  projectiles. 

1.  GUN  SYSTEM 

All  of  the  impact  tests  for  this  program  were  conducted  on  the  gun 
range  shown  in  Figure  23.  This  gun  range  is  located  in  the  Impact 
Mechanics  Facility  of  the  Air  Force  Materials  Laboratory.  For  these  tests, 
the  gun  was  used  in  two  different  configurations.  As  shown  in  Figure  23, 
the  projectile  was  explosively  driven  during  early  testing.  In  later 
tests,  the  projectile  was  driven  by  compressed  air. 

a.  Powder  Gun 

Figure  23  shows  the  original  powder  gun  configuration.  During 
launch,  the  projectile  was  housed  in  a sabot  made  of  high  density 
polyethylene.  This  heavy  sabot  was  needed  to  withstand  the  high 
acceleration  forces  generated  by  the  exploding  gas.  As  the  sabot 
traveled  down  the  gun  barrel,  the  high  pressure  gases  behind  the  sabot 
were  vented  in  the  blast  chamber.  Complete  venting  of  the  high  pressure 
gases  was  necessary  to  establish  a constant  velocity  in  the  projectile 
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Figure  23.  Overall  View  of  the  Gun  Range  Facility  Used  for  Bird 
Testing. 


prior  to  impact.  When  the  sabot  and  projectile  reached  the  end  of  the 
blast  tank,  the  sabot  was  trapped  by  the  sabot  stopper,  shown  in  Figure 
24.  The  stopper  was  so  designed  to  allow  the  projectile  to  continue  its 
flight  without,  loss  of  velocity. 

The  high  acceleration  forces  during  early  launch  and  the  rapid 
stopping  of  the  sabot  often  damaged  the  projectiles.  Requirements  for 
a thick  sabot  wall  also  limited  projectile  diameter  to  about  5 cm. 

For  these  reasons,  it  was  decided  to  redesign  the  gun  for  compressed 
air.  A more  complete  description  of  the  powder  gun  can  be  found  in  a 
report  by  Barber,  et  al.  (Reference  65).  Also  included  in  that  report 
is  a detailed  explanation  of  the  velocity  measuring  system  and  the 
photographic  system. 
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Figure  24.  Sabot  Stopper  for  Use  with  Powder  Gun. 
b.  Compressed  Air  Gun 

During  the  final  phases  of  this  project,  projectiles  were  accel- 
erated by  compressed  air  vented  from  a large  tank.  There  were  many 
advantages  to  this  system.  The  heavy,  polyethylene  sabots  were 
replaced  by  light,  thin-walled  balsa  wood  sabots.  This  was  due  to  the 
lower  acceleration  forces  Inherent  In  compressed  air  guns.  Figure  25 
Is  a picture  of  a typical  projectile  seated  in  a balsa  wood  sabot. 

In  order  to  separate  the  projectile  from  the  sabot,  a sabot  stripper  as 
shown  In  Figure  26  was  used.  In  this  technique,  the  last  section  of  the 
gun  barrel  has  a decreasing  Inner  diameter.  As  the  sabot  enters  this 
section,  it  Is  gently  decelerated  over  a distance  of  several  meters. 

This  slow  deceleration  enhances  a clear  separation  of  the  projectile 
from  the  sabot.  Barber,  et  al.  (Reference  66)  presented  the  compressed 
air  launching  system  In  greater  detail. 

2.  HOPKINSON  BAR 

In  1914,  Hopklnson  (Reference  30)  reported  the  use  of  a circular 
bar  to  measure  the  temporal  distribution  of  force  generated  during  the 
Impact  of  a lead  bullet  on  the  end  face  of  the  bar.  Previous  workers, 
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suspending  a bar  as  a pendulum,  had  measured  only  the  total  momentum 
transferred  to  the  bar  during  an  Impact  on  an  end  face.  Hopklnson 
developed  an  Ingenious  method  using  two  bars  suspending  as  pendula.  In 
which  one  bar  was  Impaated  by  the  bullet  and  the  other  acted  as  a 
momentum  trap.  Upon  Impact  of  the  first  bar,  a stress  wave,  generated  at 
the  Impact  surface,  propagated  down  the  length  of  the  bar.  This  wave 
then  passed  Into  the  second  bar  where  It  was  "trapped"  and  measured  by 
the  height  of  swing  of  the  second  bar.  By  this  method  Hopklnson  was  able 
to  resolve  the  transferred  momentum  (the  Impulse)  In  time,  thereby 
obtaining  the  force  generated  during  Impact  as  a function  of  time. 

Although  this  method  worked  well.  It  required  many  Impact  tests 
using  a like  projectile  and  like  velocity  In  order  to  resolve  the  force- 
time distribution.  In  the  present  research  programs,  the  stress  waves 
were  measured  using  strain  gages. 

a . Theory 

When  a soft  body  projectile  impacts  a Hopklnson  bar,  a shock 
wave  propagates  Into  the  projectile  as  described  earlier.  At  the  same 
time,  a stress  wave  Is  generated  in  the  bar.  It  is  essential  for  a 
Hopklnson  bar  test  that  the  pressures  generated  at  the  impact  surface 
be  much  less  than  the  yield  strength  of  the  bar  material.  For  this  case, 
the  stress  wave  in  the  bar  will  take  the  form  of  an  elastic  wave.  As  this 
wave  propagates  down  the  bar,  the  radial  stresses  will  be  released  rapidly, 
so  that  the  wave  will  become  a plane  stress  wave,  traveling  at  the 
longitudinal  wave  speed,  ca,  of  the  bar  material. 

It  Is  essential  that  the  stress  wave  in  the  bar  be  elastic  for 
several  reasons.  An  elastic  wave  will  propagate  essentially  undisturbed 
Into  the  bar.  Also,  the  strain  measured  by  the  strain  gages  will  be 
linearly  related  to  the  stress  at  that  section  of  the  bar. 

The  strain  gages  are  normally  located  at  least  ten  diameters 
down  the  bar  from  the  Impact  end.  This  Is  necessary  to  give  the  elastic 
wave  sufficient  time  to  become  planar  so  that  the  average  strain  In  the 
cross  section  will  be  the  same  as  that  measured  by  the  gages.  For  this 
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case,  the  strain  will  be  linearly  related  to  the  average  stress  at  the 
cross  section  and  the  stress  will  be  related  to  the  force  by  the  area 
of  the  cross  section. 

Figure  27  Is  an  Illustration  of  a Hcpklnson  bar.  When  the 
projectile  Impacts  the  end  of  the  bar,  the  amplitude  of  the  pressures 
generated  vary  over  the  Impact  surface.  The  resultant  stress  wave  In 
the  bar  Is  nonplanar.  However,  by  the  time  the  wave  reaches  the  strain 
gages,  It  has  become  planar  and  the  stress  measured  by  the  gages  Is  only 
the  average  of  the  Initial  stress  distribution.  Since  the  exact  stress 
distribution  cannot  be  measured,  only  the  total  force  Is  used.  The 
relation  for  force  is 

F * EAe  (62) 

where  E is  the  modulus  of  elasticity  of  the  bar,  A Is  the  cross-sectional 
area  of  the  bar,  and  e is  the  strain  recorded. 


PROJECTILE 


Figure  27,  Schematic  of  Hopkinson  Pressure  Bar. 
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Because  of  the  large  diameter  of  the  bar,  high  frequency  com- 
ponents of  the  stress  wave  are  unable  to  propagate  down  the  bar  without 
being  dispersed  (Reference  11).  Thus,  the  measured  force  profile  is 
somewhat  different  from  the  actual  forces  generated  at  the  impact  surface. 
However,  even  though  the  force  history  will  be  somewhat  distorted,  the 
total  impulse  imparted  to  the  bar  can  still  be  obtained  by  integrating 
the  force-time  trace. 


b.  AFML  Facility 

Figure  28  is  a picture  of  the  aluminum  Hopkinson  bar  mounted  in 
the  AFML  range. 


Figure  28.  Overall  View  of  the  Hopkinson  Bar  Apparatus. 


The  bar  is  7.6  cm  in  diameter  and  370  cm  long.  The  gages  are  mounted 
76  cm  from  the  impact  end.  The  bar  is  mounted  in  Teflon  rings  to  allow 
for  free  movement  and  expansion  of  the  bar.  This  is  necessary  to  prevent 
extraneous  wave  reflections  in  the  bar. 
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The  Kopklnson  bar  was  used  In  three  modes.  For  normal  Impact, 
the  end  face  of  the  bar  was  cut  at  a right  angle  to  the  axis  of  the  bar 
and  the  projectile  trajectory  was  along  the  axis  of  the  bar.  For  oblique 
Impact,  the  end  face  of  the  bar  was  cut  at  an  angle  of  25e  and  *5°  to  the 
axis  of  the  bar  and  the  projectile  trajectory  was  again  along  the  axis 
of  the  bar. 

c.  Data  Reduction 

Equation  62  was  used  to  obtain  the  force  from  the  strain  measure- 
ments. The  strain  data,  recorded  with  oscilloscopes,  was  converted  to 
force,  digitized,  and  integrated  with  time  numerically  to  obtain  the  total 
Impulse  delivered  during  impact. 

3.  PRESSURE  PLATE 

As  stated  e«rlier,  the  response  of  the  Hopkinson  bar  is  limited  by 
geometrical  wave  dispersion.  Thus,  the  force  history  recorded  using  this 
device  is  somewhat  inaccurate.  Also,  it  is  not  capable  of  measuring  the 
distribution  of  pressure  at  the  impact  surface.  The  obvious  method  for 
obtaining  this. information  is  to  mount  pressure  transducers  on  the  surface 
of  a rigid  plate.  A method  similar  to  this  was  used  by  such  researchers 
as  Bowden  and  Brunton  (Reference  67)  and  Klnslow,  et  al.  (Reference  39) 
to  measure  the  pressures  generated  during  the  impact  of  a water  jet. 
However,  as  seen  in  the  Literature  Survey,  the  results  of  previous 
researchers  using  similar  methods  of  pressure  measurement  for  the  impact 
of  water  appear  to  be  incorrect.  The  main  source  of  error  was  probably 
insufficient  frequency  response  of  the  transducers  to  measure  the  short 
duration  shock  pressures  (on  the  order  of  a microsecond  or  less).  This 
problem  was  compounded  by  the  inability  to  control  the  exact  size  and 
shape  of  the  projectile. 

Both  of  these  problems  were  overcome  in  this  study  by  using  much 
larger  projectiles  with  sufficient  material  strength  to  hold  their  shape. 
The  larger  diameter  and  the  flat  impacting  surface  of  the  projectiles 
resulted  in  shock  durations  of  10  us  or  greater  at  the  center  of  impact. 
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a.  AFML  Pressure  Plate  Facility 

Figure  29  is  a picture  of  the  pressure  plate  used  at  the  AFML 


fac  i 1 i ty . 


Figure  29.  Pressure  Plate  Apparatus. 


The  pressure  transducers  used  for  these  tests  were  piezoelectric  quartz 
pressure  transducers  manufactured  by  PCB  Corporation.  The  four  trans- 
ducers were  mounted  flush  to  the  surface  and  located  at  the  center  of 
impact  and  1.27  cm,  2.54  cm,  and  3.81  cm  from  the  center.  The  steel 
plate,  15.25  cm  in  diameter  and  5.08  cm  thick,  was  mounted  as  shown  in 
the  figure.  The  transducers  are  described  in  a report  by  Barber,  et  al. 
(Reference  65)  to  a much  greater  extent. 

b.  Data  Reduction 

Several  examples  of  actual  pressure  traces  are  shown  in  the  next 
section.  These  traces  are  digitized  from  the  original  oscilloscope 
traces.  The  amplitudes  of  the  pressures  were  obtained  by  multiplying  the 
output  voltage  of  the  gage  by  the  gage  calibration  factor. 
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The  transducers  were  not  originally  designed  for  direct  impact 
testing.  This  fact  resulted  in  several  problems  during  their  use.  The 
gages  did  not  have  adequate  acceleration  compensaticn.  Thus,  plate 
vibrations  caused  spurious  noise  to  appear  in  many  of  the  tests  results. 
Secondly,  the  resonant  frequency  of  the  gages  was  approximately  300  kHz. 

For  the  strong  shocks  which  displayed  high  amplitudes  and  short  rise  times, 
this  condition  caused  resonance  to  appear  in  the  output  of  the  gages  as 
well  as  apparent  overshoot  of  the  Initial  pressure  pulses.  This  over- 
shoot is  apparent  in  several  of  the  actual  pressure  traces  displayed  In 
the  next  section.  The  accuracy  of  the  shock  pressure  measurements  was 
also  limited  by  the  finite  frequency  response  of  the  gages  (reportedly, 

100  khz)  which  prevented  measurements  of  rise  times  of  less  than  5 us. 
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SECTION  V 

RESULTS  AND  DISCUSSION 

In  order  to  check  the  applicability  of  the  soft  body  theory  presented 
earlier,  a series  of  tests  were  conducted  using  the  Hopklnson  bar  and  the 
flat  plate  apparatus.  The  projectiles  for  the  majority  of  the  tests 
consisted  of  birds  (chickens)  and  the  bird  substitute  materials  RTV  rubber 
and  gelatin.  The  various  material  properties  used  In  the  theoretical 
analyses  are  presented  in  Appendix  C. 

All  of  the  projectiles,  except  the  birds,  were  right  circular 
cylinders.  This  greatly  simplified  comparison  of  theory  arsti  experiment. 
Both  the  RTV  and  gelatin  were  tested  with  and  without  porosity  in  order 
to  see  its  Influence  on  the  pressures  generated  during  impact.  Neoprene 
and  beef  steak  were  tested  briefly  In  order  to  study  the  effect  of  material 
strength  and  homogeneity.  The  effects  of  impact  obliquity  were  studied 
by  conducting  tests  at  three  angles,  25°,  45°,  and  90°  (normal)  to  the 
surface  of  the  plate. 

1.  HOPKINSON  BAR  TESTS 

The  assumption  that  the  various  soft  body  projectiles  behave  as 
fluids  during  impact  requires  that  the  projectile  rebound  velocity  be 
negligible.  For  normal  (90°)  impacts,  this  Implies  that  the  impulse 
imparted  to  the  target  must  be  equal  to  the  initial  momentum  of  the 
projectile,  assuming  the  target  is  rigid.  For  oblique  impact  cases, 
thi r implies  that  the  impulse  must  be  equal  to  the  normal  component  of 
the  initial  momentum,  that  is, 

I = M uo  sin  a (63) 

In  order  to  obtain  an  accurate  measure  of  the  total  impulse,  various 
soft  body  projectiles  were  impacted  against  the  Hopkinson  bar  apparatus 
described  previously.  Projectile  mass,  impact  velocity,  and  impact  angle 
were  varied  for  each  of  the  materials  tested.  The  results  of  these  tests 
are  shown  in  Figures  30,  31,  32,  and  33. 
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NORMAL  COMPONENT  OF  MOMENTUM 
mile  tin*  (N-t) 


Figure  31.  Relationship  Between  the  Normal  Component  of 

Momentum  and  the  Impulse  Imparted  to  the  Target 
for  the  Impact  of  Porous  RTV. 


NORMAL  COMPONENT  OF  MOMENTUM 
muotina  (N-t) 


Figure  32.  Relationship  Between  the  Ne  wel  Component  of  Momentum 
and  the  Impulse  Imparted  to  the  Target  for  the 
Impact  of  Por'  s Gelatin. 
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NORMAL  COMPONENT  OF  MOMENTUM, 
mu,*In«  (N-s) 

Figure  33.  Relationship  Between  the  Normal  Component  of 

Momentum  and  the  Impulse  Imparted  to  the  Target  for 
the  Impact  of  Birds. 

figure  30  presents  examples  of  the  measured  force  histories 
from  which  the  impulse  measurements  were  obtained  by  integration.  These 
are  only  approximate  force  histories  since  the  Hopklnson  bar  effectively 
filters  the  higher  frequency  components.  The  Impulse  measurements  are 
compared  with  the  normal  component  of  momentum  In  figures  31,  32,  and  33 
for  porous  RTV,  porous  (micro-balloon)  gelatin,  and  birds.  From  each 
of  these  figures.  It  is  apparent  that  there  Is  no  significant  deviation 
from  the  relationship  given  by  Equation  63  for  each  of  the  materials 
tested. 

2.  PRESSURE  PLATE  TESTS 
a.  Normal  Impact 

T 

The  fluid  model  of  soft  body  Impacts  against  a rigid  target, 
presented  earlier,  has  several  characteristic  features.  For  the  impact 
of  a cylindrical  projectile  the  pressure  history  at  the  projectile-target 
interface  Is  divided  into  three  distinct  phases;  the  shock,  the  release. 
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and  steady  flow.  For  normal  Impact,  the  Initial  Hugonlot  shock  pressure 
Is  given  by  Equation  4,  * pousuo*  anc*  ^as  a duration  at  the  center  of 

Impact  given  by  Equation  7,  tg  - a/c^,.  During  the  release  phase,  radial 
flow  Is  established.  Because  of  the  complicated  release  process,  no 
analytical  expressions  have  been  derived  for  this  phase.  However,  once 
flow  Is  completely  established,  steady  flow  relations  hold.  The  pressure 
at  the  stagnation  point,  Ps>  located  at  the  center  for  normal  Impact,  Is 
derived  from  Bernoulli's  Equation  and  takes  the  form  of  Equation  16, 


+ P 

° dP  . 

p 


2 


and  the  pressure  distribution  across  the  surface  may  be  approximated  by 
Equation  27, 


p = P {1-3  (r^r)2  + 2(rir)3) 


The  duration  of  the  entire  Impact  process  may  be  approximated  by  the 
time  that  It  takes  for  the  projectile  to  flow  through  Itself,  given  by 
Equation  21,  tQ  = L/uq. 

In  order  to  experimentally  determine  the  pressure  histories  for 
various  soft  body  impacts,  the  pressure  plate  apparatus  described 
previously  was  used.  All  of  the  projectiles,  except  the  birds,  were  right 
circular  cylinders,  with  length- to-diameter  ratios  of  approximately 
two.  During  this  testing  program,  the  size,  mass,  density,  and  porosity 
of  the  projectile  were  varied  as  well  as  the  Impact  velocity  and  Impact 
angle. 


Typical  pressure  histories  recorded  at  the  center  of  impact  for 
the  various  soft  body  materials  are  presented  In  Figure  34.  Here  the 
pressure  measured  with  the  center  transducer  (stagnation  point)  and  the 
time  are  plotted  In  the  nondlmenslonallzed  form,  obtained  by  dividing  the 
prassure  by  P = 1/2  pzuQ  (the  steady-flow  stagnation  pressure  for  an 
incompressible  fluid)  and  the  time  by  tQ  a L/uQ  (the  theoretical  duration 
of  impact).  pz,  the  average  density  of  the  porous  material  was  deflnad 
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In  Equation  58.  The  impact  velocity  for  each  of  these  tests  was 
approximately  200  m/s.  The  similarities  of  the  pressure  histories  for 
the  various  materials  are  evident.  Each  shows  the  expected  fluid 
characteristics;  namely,  an  initial  peak  pressure  followed  by  a decay  to 
a lower  flow  pressure,  with  a total  duration  of  approximately  tQ  = L/uQ. 

Of  course,  the  exact  shape  of  the  pressure  histories  varies  for  each 
material  due  to  differences  in  material  properties. 

The  peak  pressure  is  seen  to  be  greatest  for  the  pure  gelatin, 
due  to  its  high  shock  velocity  and  lack  of  porjsity.  Each  of  the  other 
materials,  with  the  exception  of  neoprene,  had  certain  amounts  of  porosity 
which  greatly  decreased  their  shock  velocities,  hence  their  shock  pres- 
sures. Neoprene,  even  with  zero  porosity,  has  a very  low  shock  velocity. 

On  some  of  the  bird  shots  the  pressure  did  not  rise  immediately  to  a peak. 
This  was  due  to  the  impact  of  such  things  as  the  bird's  feet  or  feathers 
prior  to  impact  by  its  main  body. 

For  each  of  these  materials,  an  apparent  steady  flow  region 
exists  and  it  is  most  readily  seen  in  the  pressure  history  of  porous  RTV. 
Each  of  the  other  materials  displayed  a tendency  for  the  pressure  to 
slowly  decay  in  this  flow  regime.  This  decay  from  steady  flow,  most 
obvious  for  the  beef  and  pure  gelatin,  was  described  in  the  theory  as 
being  the  unsteady  flow  effects  due  to  the  finite  length  of  the  projectile. 
As  stated  in  the  theoretical  section,  this  effect  should  be  greatest  for 
subsonic  impacts  such  as  those  of  the  beef  and  gelatin.  It  should  be 
least  apparent  for  supersonic  impacts  such  as  for  the  porous  RTV  and 
micro-balloon  (porous)  gelatin. 

Since  pure  gelatin  Is  almost  incompressible  at  this  impact  velo- 

p 

city,  its  stagnation  pressure  is  'een  to  be  approximately  1/2  pzuq  . 

Each  of  the  other  materials  showed  varying  degrees  of  compressibility, 
evidenced  by  the  fact  that  in  the  steady  flow  regime,  their  nondimension- 
all  zed  pressures  had  values  greater  than  one. 

One  of  the  major  features  distinguishing  the  pressure  histories 
of  birds  from  that  of  porous  RTV  and  micro-balloon  gelatin  is  the  much 
larger  amount  of  high-frequency  "noise"  superimposed  on  the  pressure 
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profile  of  birds.  Initial  inquiries  into  the  source  of  this  noise  were 
fruitless.  In  order  to  see  if  it  was  caused  by  the  bone  structure  of  the 
bird,  a beef  projectile  was  constructed  and  tested.  As  can  be  seen, 
the  noise  was  even  greater  for  the  beef  than  for  the  birds.  Thus,  the 
effect  of  the  bone  structure  was  assumed  negligible.  Another  thought 
was  that  the  noise  might  have  been  caused  by  the  break-up  or  tearing  of 
the  material  (creation  of  new  surfaces)  during  impact.  This  theory 
seemed  resonable  since  both  the  birds  and  the  beef  consisted  of  muscle 
fibers  which  would  require  more  energy  to  tear  than  would  the  RTV  or 
gelatin.  Along  this  line  of  thought,  neoprene  was  tested  since  it  bounces 
at  velocities  of  200  m/s  and  there  is  no  resulting  material  break-up. 
Accordingly,  the  pressure  history  of  the  neoprene  impact  showed  no  noise 
except  that  due  to  acceleration  loads  on  the  gages  caused  by  the  pressure 
plate  vibrations.  Thus,  there  appears  to  be  a direct  correlation  between 
the  amount  of  energy  expended  in  material  break-up  and  the  amount  of 
high-frequency  noise  recorded.  (The  extreme  amount  of  noise  recorded 
during  the  impact  of  pure  gelatin  was  apparently  due  to  excitation  of  the 
gage  near  its  resonant  frequency). 

A point  of  interest  noted  during  this  testing  program  was  the 
fact  that  the  neoprene's  pressure  history  was  somewhat  similar  to  that  of 
the  other  materials  even  though  it  did  not  "flow"  at  impact  velocities 
as  high  as  200  m/s.  This  was  noted  for  porous  RTV  and  micro-balloon 
gelatin  which  bounced  at  impact  velocities  of  100  m/s.  This  would  seem 
to  imply  that  the  basic  fluid  flow  theory  presented  in  this  report  holds 
even  for  materials  which  do  not  flow  (or  shatter)  during  impact  but 
which  do  undergo  large  amount  of  deformation  (strain). 

A closer  comparison  of  the  pressure  histories  for  the  impact  of 
birds  and  the  two  bird  substitute  materials,  porous  RTV  and  micro-balloon 
gelatin,  can  be  obtained  from  Figures  35,  36,  and  37.  In  these  figures, 
typical  pressure  histories  at  the  center  of  impact  are  shown  for  various 
impact  velocities.  In  each  figure  a short  line  representing  the  pre- 
dicted shock  pressure  amplitude,  PH>  and  duration,  tg,  is  shown  as  well 
as  a long  line  representing  the  predicted  amplitude  of  the  steady  flow 
pressure  at  the  stagnation  point,  P . The  pressure  and  time  are  again 
presented  in  nondimensional ized  form. 
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Figure  37.  Typical  Pressure  Data  at  Various  Impact  Velocities 
for  Birds,  Normal  Impact. 
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The  calculated  amplitude  and  duration  of  the  shock  pressures 
agree  well  with  the  measurements  for  the  two  bird  substitute  materials, 
although  some  transducer  resonant  excitation  and  resultant  overshoot  Is 
seen  on  several  shots.  Discrepancies  between  theory  and  measurements 
for  birds  should  be  expected  due  to  the  rounded  qeometry  and  Irregular 
surface  of  the  bird  and  the  prior  Impact  of  su.  *-hings  as  feathers, 
wings,  and  legs. 

As  was  stated  earlier,  each  of  the  materials  seems  to  display  a 
regime  of  steady  flow  during  which  the  measured  pressure  agrees  with  the 
calculated  stagnation  pressure.  The  deviations  in  the  pressure  histories 
for  the  micro-balloon  gelatin  at  high  velocity  may  have  been  due  to 
geometric  distortion  of  the  projectile  during  launch  or  to  resonant 
excitation. 

In  order  to  demonstrate  the  radial  variation  of  pressure, 

Figure  38  presents  a typical  porous  RTV  Impact  in  which  the  pressure 
histories  at  three  transducer  locations  are  shown.  The  peak  pressure 
appears  to  be  less  at  the  transducer  located  at  0.6  r/a  than  at  the 
center  since  the  duration  of  the  shock  pressure  at  this  point  Is  so  short 
that  the  transducer  could  not  respond  rapidly  enough  to  record  it.  Of 
course,  there  was  no  shock  at  the  transducer  located  outside  the  initial 
impact  area  of  the  projectile.  This  figure  also  demonstrates  the  radial 
decrease  in  pressure  In  the  steady  flow  regime. 

Figure  39  shows  a typical  pressure  history  for  porous  RTV 
recorded  along  the  major  axis  of  an  oblique  Impact  of  45°.  Lines  are 
drawn  in  each  figure  representing  the  calculated  shock  pressure,  and 
full,  steady-flow  stagnation  pressure.  The  predicted  shock  durations  are 
not  shown  In  these  figures.  As  explained  in  the  theory,  the  shock 
pressure  for  an  Impact  velocity  uQ  and  an  impact  angle  a should  be 
identical  to  that  produced  by  a normal  Impact  with  an  initial  velocity 
of  uQ  sin  a.  From  Equation  36,  the  45°  impact  angle  Is  seen  to  be 
very  near  the  critical  angle  for  RTV.  This  Implies  that  the  full  shock 
pressure  may  exist  over  the  entire  Impact  surface  and  that  the  duration 
of  the  shock  near  the  center  of  impact  should  be  sufficient  for  the 
pressure  gages  to  record  the  full  shock  pressure.  Also,  expected  from 
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Figure  38.  Typical  Distribution  of  Pressure  PUte  Data  for  the 
Impact  of  Porous  RTV,  Normal  Impact. 
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Figure  39.  Typical  Distribution  of  Pressure  Plate  Data  for  the 
Impact  of  Porour  RTV,  Oblique  Impact  ( a ■ 45*). 
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the  theory  is  the  fact  that  the  stagnation  point,  which  is  approximately 
3 cm  above  the  center  transducer,  should  "see"  the  full  stagnation  pressure 
during  steady-flow,  regardless  of  the  angle  of  impact. 

From  the  figure,  several  of  these  points  are  clear.  The  shock 
pressure  at  the  center  gage  Is  the  full  shock  pressure,  and  the  other 
gages  record  peak  pressures  near  this  value.  The  steady-'flow  pressure  is 
also  seen  to  be  highest  at  the  gage  2.54  cm  above  center.  The  fact  that 
the  pressure  here  is  somewhat  below  the  stagnation  pressure  implies  that 
this  gage  was  not  located  exactly  at  the  stagnation  point.  It  should  be 
noted  that  for  much  smaller  angles  of  impact  (greater  impact  obliquity' 
the  amplitude  of  the  calculated  shock  pressure  will  be  less  than  the 
steady-flow  stagnation  pressure  for  RTV. 

In  Figures  34-39,  typical  pressure  histories  for  the  various 
materials  were  presented  in  order  to  demonstrate  a qualitative  under- 
standing of  the  pressures  generated  during  soft  body  impact.  However, 
in  order  to  compare  the  fluid  theory  with  the  soft  body  impact  experi- 
ments, a more  quantitative  presentation  of  the  pressure  data  is  required. 

In  this  section,  the  experimental  results  are  compared  with  theoretical 
predictions  for  the  various  materials. 

The  variation  of  shock  pressure  with  impact  velocity  for  normal 
impact  is  given  in  Figures  40,  41,  and  42  for  porous  RTV,  micro-balloon 
gelatin,  and  birds,  respectively.  In  Figure  40,  the  theoretical 
Hugoniot  shock  pressure  is  given  for  both  pure  (z  = 0.00)  and  porous 
(z  = 0.50)  RTV.  Experimental  measurements  for  the  porous  RTV  tend  to 
be  higher  than  that  predicted.  However,  as  stated  previously,  it 
appears  in  Figure  36  that  the  high  peak  pressures  may  have  been  due  to 
some  transducer  overshoot.  In  fact,  the  pressures  tended  to  oscillate 
about  the  predicted  shock  pressures.  In  Figure  41,  the  theoretical 
Hugoniot  shock  pressure  is  given  for  pure  (z  = 0.00)  and  porous  (z  * 0.40) 
gelatin.  Experimental  measurements  for  the  Impact  of  the  micro-balloon 
gelatin  are  shown  to  agree  very  closely  with  the  predictions. 

In  Figure  42,  the  theoretical  Hugoniot  shock  pressure  is  given 
for  birds.  As  stated  in  Appendix  C,  this  curve  is  based  on  the  assumption 
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Figure  41.  Hugonlot  Shock  Pressures 
of  Porous  Gelatin.  (1  N 
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Figure  42.  Hugonlot  Shock  Pressures  Measured  During  Normal  Impact  of 
Birds.  (1  MN/m2  - 10  Bars,  103  kg/m3  ■ .036  lb/1n3). 
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that  the  Hugonlot  relationship  for  birds  is  similar  to  that  for  water 
(or  gelatin)  with  an  assumed  pure  density  of  * 1060  kg/m  and  a 
porosity  of  z * 0.10  (average  density  p2  ■ 950  kg/m3).  The  data  shown 
here  were  obtained  from  tests  conducted  on  small  birds  at  AFML  and  on 
large  birds  at  AEDC  (Reference  66).  The  measured  peak  pressures  fall 
below  that  predicted.  The  fact  that  the  peak  pressure  tended  to  Increase 
with  bird  mass,  hence  with  bird  size,  Implies  that  the  low  peak  pressure 
measurements  may  be  due  to  the  lack  of  planar  Impact  for  birds.  As  shown 
In  the  theoretical  section,  the  duration  of  the  full  shock  pressure  Is  a 
function  of  the  radius  of  curvature  of  the  Impacting  surface,  so  that  the 
shock  duration  at  the  center  of  Impact  Is  much  less  for  a projectile  with 
a rounded  end  than  for  a right  circular  cylinder.  For  small  birds  with 
an  Irregular  impact  surface  and  a small  radius  of  curvature,  the  shock 
duration  was  probably  so  small  that  the  transducers  could  not  respond 
rapidly  enough  to  record  the  shock  pressure.  For  larger  birds  with  larger 
radii  of  curvature,  the  duration  of  the  shock  was  greater.  The  peak 
pressures  measured  on  several  of  the  4 kg  bird  tests  approached  the 
predicted  values.  Other  reasons  for  the  low  peak  pressures  Include  the 
existence  of  material  such  as  feathers  and  legs  which  Impacted  prior  to 
the  main  body  and  the  existence  of  lower  impedance  material  at  the 
surface  of  the  main  body. 

The  variation  of  pressures  in  the  steady  flow  regime  with  impact 
velocity  Is  presented  In  Figures  43,  44,  and  45.  The  pressure  at  the 
center  stagnation  point  is  plotted  for  RTV  in  Figure  43.  Theoretical 
predictions  are  given  for  pure  and  porous  RTV  using  the  compressible 

p 

fluid  theory.  A curve  of  P a 1/2  p_u_  (Incompressible  fluid  theory, 

3 1 0 

P2  = 670  kg/m  ) Is  also  presented  for  comparison.  The  measured  stagnation 
pressures  for  porous  RTV  agree  very  closely  with  those  predicted  by 
the  compressible  theory. 

In  Figure  44,  the  steady-flow,  stagnation  pressure  is  given  for 
both  pure  gelatin  and  micro-balloon  gelatin  impact  tests.  Theoretical 
curves  *rom  compressible  theory  arc  shown  as  well  as  the  plot  of 
P = 1/2  pzuq  (p2  = 640  kg/m'  ).  The  experimental  data  is  seen  to  agree 
very  nicely  with  the  compressible  theory  for  both  types  of  gelatin. 
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Figure  43.  Steady  Flow  Pressures  Measured  at  the  Center  of  Impact 
During  Normal  Impact  of  Porous  RTV.  (1  MN/nr  ■ 10  Bars, 
103  kg/m*  - .036  lb/1n3). 
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Figure  44.  Steady  Flow  Pressures  Measured  at  the  Center  of  Impact 

During  Normal  Impact  of  Porous  Gelatin.  .(1  MN/nr  * 10  Bars, 
103  kg/m3  ■ .036  lb/1n3). 
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Figure  45.  Steady  Flow  Pressures  Measured  at  the  Center  of  Impact 
During  Normal  Impact  of  Birds.  (1  NN/m2  » 10  Bars, 

103  kg/m3  ■ .036  lb/1n3). 
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In  Figure  45,  the  stagnation  pressure  Is  given  for  the  Impact  of 
birds.  A theoretical  prediction  Is  made  assuming  Pf  * 1060  kg/m  and 
z ■ 0.10  and  ■ 950  kg/m^  and  z ■ 0.0.  As  shown  earlier  In  Figure  35 
for  the  Impact  of  birds,  end  effects  In  the  bird  causes  a gradual  decrease 
In  flow  pressure,  making  measurement  of  the  steady  flow  pressure  difficult. 
However,  the  measured  data  does  appear  to  be  In  fair  agreement  with  the 
theoretical  curve  based  on  the  compressible  flow  theory,  although  there 
Is  more  scatter  In  the  data  than  seen  previously  for  the  RTV  and  gelatin. 

Figures  46,  47,  48,  and  49  show  the  variation  of  pressure  In  the 
steady  flow  regime  with  radial  distance  from  the  center.  The  pressure  Is 
again  nondlmenslonallzed  and  the  radial  distance  Is  nondl mens Iona 11 zed 
with  respect  to  the  Initial  radius  of  the  projectile,  a. 
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Figure  46.  Radial  Distribution  of  Steady  Flow  Pressures  for  the 
Normal  Impact  of  Porous  RTV. 
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(ADJUSTED) 

Figure  49.  Radial  Distribution  of  Steady  Flow  Pressures  for  the 
Normal  Impact  of  Birds. 

The  theoretical  curves  are  drawn  from  the  approximate  relationship  given 
by  Equation  27.  The  experimental  data  are  presented  In  the  form  of  bars 
which  designate  the  data  bounds.  The  large  scatter  In  the  data  may  have 
been  due  to  the  projectiles  hitting  off-center  and  Inaccuracies  In  the 
pressure  measurement.  Because  of  the  data  scatter.  It  Is  not  possible  to 
determine  how  accurately  the  theoretical  curves  predict  the  actual  pressure 
distribution. 

b.  Oblique  Impact 

Original  plans  called  for  experimental  testing  of  normal  Impact 
only.  However,  during  the  course  of  this  research  program.  It  became 
evident  that  some  Investigation  of  oblique  Impact  should  be  performed 
In  order  to  test  a few  of  the  basic  concepts. 

As  already  noted  In  the  Hopklnson  bar  tests,  It  was  shown  that 
the  Impulse  Imparted  to  the  target  was  related  to  the  Initial  momentum  of 
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the  projectile  by  the  sine  of  the  Impact  angle;  thet  Is,  the  Impulse  was 
•quel  to  the  norms!  component  of  momentum.  As  wes  demonstrated  In 
Figure  39,  the  aaqtlltude  of  the  Hugonlot  shock  pressure  for  en  oblique 
Impact  of  engle  a It  simply  the  shock  pressure  predicted  for  the  normel 
component,  uQ  sin  a,  of  the  Impect  velocity.  For  steed?  flow,  the 
stegnetlon  pressure  Is  the  seme  es  for  normel  Impect,  elthough  the 
stegnetl on  point  Is  no  longer  loceted  et  the  center  of  Impect,  end  the 
distribution  of  pressure  over  the  surfece  Is  no  longer  exl symmetric. 

For  the  purpose  of  this  limited  testing,  only  the  Impect  of 
porous  RTV  wes  considered.  Impect  tests  were  conducted  et  engles  of  25* 
end  45*  to  the  surfece  of  the  terget  plete.  In  Figure  39,  typlcel  pres- 
sure histories  elong  the  mejor  axis  for  en  oblique  Impect  of  45*  were 
given.  In  Figures  50  end  51 , the  celculeted  Hugonlot  shock  pressures 
ere  compered  with  the  meesured  peek  pressures  for  oblique  Impects  of  45* 
end  25°,  respectively.  The  pressure  et  eech  of  the  four  trensducers  Is 
presented  for  the  45*  Impect  case.  For  the  25*  Impect  cese,  only  the 
date  from  the  gage  nearest  the  Initial  Impact  point  Is  given.  For  this 
cese,  the  other  transducers  do  not  see  the  full  shock  pressure  since  the 
critical  engle  has  been  exceeded.  The  general  trend  of  the  data  seems 
to  agree  with  the  theory.  The  few  high  measurements  may  have  been  due 
to  electronic  overshoot  of  the  center  transducer.  This  belief  Is 
strengthened  by  the  feet  thet  this  discrepancy  only  occurred  et  the  center 
transducer  for  the  highest  velocity  Impects. 

In  Figure  52,  the  steady  flow  pressure  measured  at  the  transducer 
nearest  the  stagnation  point  Is  presented  as  a function  of  velocity. 

Date  for  both  90*  end  45*  Impect  of  porous  RTV  are  compared  with  the 
stegnetlon  pressure  predicted  by  Equation  27.  The  measurements  for  90* 
Impect  (presented  previously  In  Figure  43)  compere  very  favorably  with 
the  predictions.  However,  the  date  for  45*  Impect  Is  less  t^en  thet 
predicted.  This  was  probably  due  to  the  fact  that  the  transducer  was 
not  loceted  exactly  at  the  stagnation  point  or  thet  the  region  of  full 
stegnetlon  pressure  wes  smeller  then  the  active  area  of  the  gage. 
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Figure  50.  Hugonlot  Shock  Pressures  Measured  at  Various  Points  Alona 
the  Target  Surface  During  the  Oblique  Impact  (o  ■ 45*)  of 
Porous  RTV.  (1  MN/«r  ■ 10  Bars). 
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Figure  51.  Hugonlot  Shock  Pressures  Measured  at  the  Initial  Point  of 
Impact  During  the  Oblique  Impact  (a  ■ 25°)  of  Porous  RTV. 
(1  MN/m2  - 10  Bars). 
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Figure  52.  Steady  Flow  Pressures  Measured  at  or  Near  Stagnation  Point  for 
Normal  and  Oblique  (a  ■ 45°)  Impact  of  Porous  RTV. 

(1  MN/b>2  . io  Bars) 
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SECTION  VI 
CONCLUSIONS 

Based  on  the  experimental  results,  the  following  conclusions  were 
drawn. 


1.  An  unsteady,  hydrodynamic  model  has  been  developed  which 
describes  the  Impact  behavior  of  low  strength  (soft  body)  projectiles 
against  a flat  target. 

2.  This  model  has  been  verified  experimentally  for  normal  and 
oblique  Impacts  of  porous  and  nonporous  projectile  materials  at  Impact 
velocities  sufficient  to  cause  large  deformations  or  fluid  like  flow 
of  the  projectile. 

3.  Based  on  this  model  and  Its  subsequent  verification,  a bird 
(chickens)  can  best  be  described  as  a low  strength  material  with  the 
equation  of  state  of  water,  an  average  density  of  950  kg/m3  (.034  lbm/1n.3) 
and  10X  porosity. 
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APPENDIX  A 

RADIAL  PRESSURE  DISTRIBUTION 

The  pressure  distribution  along  the  Impact  surfaca  during  steady 
flow  may  be  approximated  by  Equation  26 


P « P,  exp  (-^(J)2) 


(A-l) 


or  by  Equation  27 


P * P_  ♦ 2(~— ) *} 


<2* 


(A-2) 


These  two  expressions  were  developed  by  Banks  and  Chandrasekhara 
(Reference  56)  and  Leach  and  Walker  (Reference  58)  respectively.  Both 
relations  must  satisfy  momentum  considerations*  given  here  by  Equation  22 

F « pAuQ2  (A-3) 

The  force.  F.  at  the  surface  Is  simply  the  Integrated  pressure,  such  that 

F * /PdA  (A-4) 

For  a cylindrical  projectile.  Equations  A-3  and  A-4  may  be  combined  to 
obtain 


2 2 /" 
np»  u ■ 2*1  P r dr  , 

0 Jo 


(A-5) 


The  constant  may  be  obtained  by  substituting  Equation  A-l  Into 
Equation  A-5  such  that 

n -‘i#’  \ 

•\  * 2 V*  <7  Jr 


P«4 


(A-6) 


Solving  this  expression  gives 


h m pu 


(A-7) 
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Tha  constant  may  bo  obtalnad  by  substituting  Eou&tlor,  A-2  iotu 
Equation  A-5 

»V  • * \ /o'J*(l-s<!7)J  * <*-•> 

Solving  this  axprasslon  glvas 

2/P.  (A-9) 

2 O C 


i 
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APPENDIX  B 

GENERAL  PRESSURE  - DENSITY  RELATIONSHIPS 


In  this  section,  the  simple  mixture  theory  given  by  Torvlk 
(Reference  64)  Is  used  to  develop  tho  pressure-density  relationships  for 
a porous  material  under  both  Isentroplc  and  shock  compression.  Tory  Ik's 
thaory  predlcttd  tha  pressure-density  relationship  across  a shock  for  a 
homogeneous  mixture  based  on  the  Hugonlot  relationships  of  the  constituents. 
In  this  theory  several  fundamental  assumptions  are  made  as  follows: 

1.  The  mixture  Is  maeroscoplcally  homogeneous  and  Isotropic. 

2.  The  density  of  the  mixture  Is  the  sum  of  the  products  of  the 
volume  fraction  and  density  of  each  component. 

3.  The  density  of  each  component  Is  the  same  as  that  predicted 
for  a homogeneous  sample  If  the  Individual  material  at  the  same 
pressure. 

4.  The  size  of  Individual  particles  of  each  component  Is  such  that 
the  shock  pressure  Is  the  same  In  each  component. 

5.  The  constituents  do  not  react  or  change  phase  during  the 
compression  process. 

From  the  second  assumption  It  Is  clear  that 


PAVG 


N 

£ 


fvi 


i=I  pi 


(B  1) 


where  Is  the  average  density  of  the  mixture,  fv^  Is  the  volume 
fraction  of  the  1th  component,  Is  the  density  of  the  1th  component, 
and  N Is  the  number  of  components.  Rewriting  Equation  8-1  In  terms  of 
the  mass  fraction 


i 

PAV6 


N 

l 

i=i 


(B-2) 


where  fm,|  Is  the  mass  fraction  of  the  1th  component. 
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Assuming  that  the  mass  fraction  of  tha  components  Is  unchanged 
aftar  tha  material  compression,  it  <s  taan  that 


H 

P1AVG  m 1«1  *L 
°2AVG  **1 


(B-3) 


where  tha  subscript  1 rafars  to  tha  Initial  (unstrassad)  stata  and  2 tha 
final  (strassad)  stata.  Equation  B-3  nay  ba  simplified  by  realizing  that 


and 

pi 

y * fv^  (k  * 1*  2,  3,...H) 

r — i- 

i»l  plj 

Thus.  Equation  B-3  can  ba  rewritten 

P1AVG  N Pli 

= E fv.  (-— ) 

P2AVG  i=l  1 P2j 


(B-4) 


(B-S) 


(B-6) 


Although  Equation  B-6  was  darlvad  by  Torvlk  for  danslty  changas  across  a 
shock.  It  should  hold  trua  for  any  compression  process  In  which  all  tha 
assumptions  art  satlsflad  and  tha  prassura-danslty  relationships  are 
known  for  each  constituent  of  the  mixture.  Therefore,  Equation  B-6 
will  ba  used  to  derive  the  pressure-density  relationships  for  both 
shock  and  Isentroplc  compression. 
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1 . ISENTROPIC 

The  pressure-density  relationship  of  tho  toft  body  material  will  bo 
represented  by  Equation  S3 

p9  « A 

2 pi 

Assuming  tho  mttoriol  con  bo  roprosontod  by  a llnotr  Hugonlot  (Equation  48) 

us  * Co  + 

tho  constants  can  be  expressed  by  Equation  54 

A = p.C  2/(4k-l) 

JL  O 

B = 4K-1 


Equation  53  may  bo  rearranged  to  obtain 

HMl/B 


(B-7) 


The  pressure-density  relationship  of  air  for  an  Isentroplc 
compression  may  bo  approximated  by  tho  relationship  for  a perfect  gas 


nr 


where  y Is  the  ratio  of  specific  heats  (y  ■ 1.4  for  air).  Thus,  sub- 
stituting Equations  B-7  and  A-6  Into  Equation  B-6,  the  pressure-density 
relationship  under  Isentroplc  compression  for  a porous  material  Is  seen 
to  be 


P1AVG 

P2AVG 


(B-9) 


where  subscript  m represents  the  soft  body  material  and  n represents  the 
air. 
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2.  SHOCK 

As  stated  In  the  section  on  equations  of  state,  for  pressures  In 
the  rang*  of  Interest,  tho  pressure-density  relationships  under  Ison* 
tropic  and  shock  compression  art  approximately  equal.  Thus,  Equation  B-7 
will  be  used  to  represent  the  shock  pressure-density  relationship  for 
soft  body  materials. 


For  air  the  shock  relation  Is  significantly  different  from  the 
Isentroplc  relation.  The  shock  relation  for  a perfect  gas  Is 


l + 


Y-l 


V+l  ^2 
Y 1 p 


(B-IO) 


This  equation  Is  plotted  In  Figure  53  along  with  experimental  values  for 
real  air  obtained  by  Deal  (Reference  68).  As  can  be  seen  from  the 
figure,  the  density  ratio  for  a perfect  gas  approaches  a fixed  value 
(6.0)  as  the  pressure  Increases.  This  Is  not  true  with  the  experimental 
data  for  air. 

In  order  to  develop  a more  realistic  equation  of  state  for  air, 
the  data  of  Deal's  work  was  studied  In  more  detail.  Figure  54  Is  a 
plot  of  the  shock  velocity  versus  particle  velocity  for  air.  From  this 
data.  It  can  be  seen  that  air  can  be  represented  by  the  linear  Hugoniot 
relationship 


t 


where  k * 1.03.  In  the  section  on  Equations  of  State  It  was  noted  that 
materials  which  could  be  represented  by  the  linear  Hugoniot  relation 
could  also  be  represented  by  the  pressure-density  relationships  given 
In  Equations  52  or  53,  Upon  examination,  air  was  best  represented  by 
Equation  52  which  has  the  form 

P « pc  2q/(l-kq)a  (B-ll) 

2 1 O 
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A plot  of  Equation  Ml  In  Flgura  S3  was  toon  to  fit  tho  exparlmantal  dot* 
vary  woll. 

Rearranging  Equation  B-ll  Into  • form  compatible  with  Equation  B-6 


f • =*  -j-i-2-  — 3 — ] 

P1  Po  (l-kq)^ 


(B-12) 


which  glvos 


P (1-kq)  - 


(B-13) 


Solving  by  tho  quadratic  aquation.  It  Is  soon  that 


whara 


. „ j 

P,0 

2?k  + ~y- 


(B-14) 


- {(2Pk  + ^~)2  - 4P2k2}y2Hk2 


IB-15) 


It  can  be  noted  from  tha  definition  of  q that 


.2  „ JL 
p.  i-q 


(B-16) 


P 2 

In  ordtr  for  r*  to  bt  posltlvt*  q must  bt  loss  thtn  unity*  For  this  to 
pi 

be  true.  Equation  B-14  must  taka  tha  form 


* ■ - qa 


(B-17 ) 


no 
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Thus*  tha  final  expression  for  the  shock  pressure-density  relationship 
of  sir  Is  given  by  Equation  B-16  where  q Is  given  by  Equation  B-17  end 
q^  end  q{  by  Equation  B-15. 

i : > ••  . . , ••  . 

Substitution  of  Equstlons  B-7  and  b-16  Into  Equation  A-4  gives  the 
pressure-density  relationship  under  shock  compression  for  a porous 
material  of  the  form 

-1^ 

« fv/-l  ^ + fv  (l-q)  (B-it) 

P2AVG  A / 

where  subscript  m represents  the  toft  body  material  and  n represents 
real  air,  and  q Is  defined  by  Equations  B-17  and  B-15. 
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APPENDIX  C 
MATERIAL  PROPERTIES 

The  properties  of  the  various  materials  studied  In  this  work  are 
presented  In  this  section. 

1 . WATER 

For  water,  Hermann  (Reference  38)  showed  that  the  linear  Hugonlot 
takes  the  form 

u = c + 2.0  u 
so  p 

where  cQ  Is  the  normal  sound  speed  In  water.  This  relation  was  shown 
to  provide  a quite  accurate  fit  to  the  exact  data  for  water  up  to 
Mach  >1.2.  The  density  and  wave  speed  for  water  at  a temperature  of 
20°C  and  pressure  of  1 atmosphere  (Reference  62)  are 

pf  = 998.23  kg/m3  (.0361  lb/in3) 

c = 1482.9  m/s  (4865  ft/a) 
o 


2.  GELATIN 

For  this  program,  commercial  gelatin  with  a strength  of  250  Bloom 
was  tested.  Historically,  gelatin  has  been  used  to  simulate  flesh  in 
bullet  Impact  tests.  For  this  reason,  and  because  of  Its  low  cost  and 
ease  of  fabrication,  this  material  is  currently  used  by  several  aircraft 
engine  blade  manufacturers  to  simulate  birds  in  blade  Impact  tests. 

During  this  test  program,  it  exhibited  sufficient  shear  strength  to 
withstand  the  high  acceleration  loads  attained  during  projectile  launch 
In  the  gun  range.  However,  difficulty  In  launch  did  occur  above  velo- 
cities of  250  m/s.  Porosity  was  attained  by  the  addition  of  phenolic 
micro-balloons. 

because  of  Its  unlikely  Involvement  In  Impact,  apparently  no 
compression  data  has  been  generated  for  gelatin.  However,  the  assumption 
made  by  previous  researchers  has  been  to  approximate  gelatin  with  the 
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•qua t Ion  of  state  of  water.  The  only  major  drawback  to  this  Is  tht  fact 
that  th«  normal  aquation  of  stata  of  water  doas  not  Incorpora to  shaar 
strangth.  Tha  llnaar  Hugonlot  Is  taken  to  ba  that  of  watar 

u9  ■ °0  + 2<0up 

with 

pf  • 1060  kg/*3  (.0383  lb/in3> 

e “ 1482.9  m/s  (4865  ft/a) 

Q 

The  micro-balloon  gelatin  was  treated  as  simply  a mixture  of  gelatin  and 
air.  For  these  tasts,  the  material  had  tha  properties 

s = .40 

PE  » 640  kg/*2  ( .0231,  lb,/ in3) 

3.  RTV  RUBBER 

RTV  (Room  Temperature  Vulcanized)  rubber  Is  a silicons  elastomer 
which  Is  presently  used  by  several  organizations  as  a bird  substitute 
material.  In  previous  work  conducted  at  the  AFML  Impact  Facility*  It 
was  determined  that  the  linear  Hugonlot  for  RTV-11  could  be  approximated 
by  the  relationship 

u ■ c ♦ 3*62  u 
no  p 

For  the  current  study,  a mixture  of  RTV-560  and  RTV-921  was  tested. 
Porosity  was  obtained  by  the  use  of  a blowing  agent.  Since  no  Hugonlot 
date  existed  for  this  material,  the  relationship  for  RTV-11  was  used. 

It  was  understood  that  use  of  this  relationship  might  cause  some 
Inaccuracies  In  the  resulting  calculations.  The  original  560/921 
material  had  a density  and  wave  speed  of 

pf  * 1,330  kg/*3  (.0481  lb/in3) 

o » 830  m/a 


(2743  ft/a) 
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With  the  blowing  agent  added,  the  porous  RTV  had  the  properties 


B a .50 

0 = 670  k«/m3  (.0242  lb/in3) 

s 

4.  BIRDS 

All  bird  tests  conducted  In  this  program  used  chickens  as  the 
projectiles.  However,  many  researchers  feel  that  chickens  are  not 
representative  of  flying  birds.  The  major  point  of  contention  Is  that 
flying  birds  have  much  lower  densities  than  chickens.  If  this  Is  the 
case,  then  the  equation  of  state  derived  In  this  section  for  chickens 
will  not  be  adequate  for  other  birds,  although  the  only  difference  will 
be  the  amount  or  assumed  porosity. 

A large  number  of  tests  In  which  baby  chickens  have  been  Impacted 
against  a rigid  plate  have  been  conducted  In  the  Impact  Mechanics 
Facility  over  the  past  several  years.  This  work  was  done  In  conjunction 
with  FOD  (Foreign  Object  Damage)  studies  of  aircraft  engine  blading 
and  aircraft  windshields.  Similar  testing  with  large  chickens  has  been 
conducted  at  the  Arnold  Engineering  Development  Center.  During  these 
studies,  the  birds  have  been  accelerated  to  velocities  exceeding  300  m/s. 

Although  birds  are  Inhomogeneous,  anisotropic,  and  not  axl symmetric 
it  was  felt  that  to  some  degree,  they  could  be  approximated  by  a cylinder 
of  material  which  was  both  homogeneous  and  Isotropic.  Thus,  the 
assumption  made  in  this  program  Is  that  they  may  be  approximated  by 
a cylinder  of  porous  water.  Here  the  shear  sti'ength  of  the  bird  has 
been  neglected. 

Thus,  the  linear  Hugonlot  is  taken  to  be  that  of  water 

u = c + 2.0  u 
so  p 

with 
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The  density  measured  for  several  chickens  of  various  sizes  was  found 
to  be 

O » 950  kf/»3  (.0342  lb/in3) 

Several  researchers*  Including  Tudor*  et  al.  (Reference  6)  have  measured 
the  density  of  various  parts  of  the  chicken  to  be  approximately 
p • 1*060  kg/m3.  If  this  Is  actually  the  case*  then  chickens  can  be 
considered  to  have  a full  density  of  ■ 1*060  kg/m3  with  10X  porosity 
(p2  • 950  kg/m3).  However*  many  researchers  feel  that  chickens  can 
be  better  approximated  by  a material  with  an  original  density  of 
Pf  ■ 950  kg/m3  with  no  porosity*  Both  of  these  models  will  be  checked 
against  experimental  results. 
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